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Abstract. Magnetars are the strongest magnets in the present universe and the 
combination of extreme magnetic field, gravity and density makes them unique 
laboratories to probe current physical theories (from quantum electrodynamics to 
general relativity) in the strong field limit. Magnetars are observed as peculiar, burst- 
active X-ray pulsars, the Anomalous X-ray Pulsars (AXPs) and the Soft Gamma 
Repeaters (SGRs); the latter emitted also three “giant flares”, extremely powerful 
events during which luminosities can reach up to 10^^ erg/s for about one second. The 
last five years have witnessed an explosion in magnetar research which has led, among 
other things, to the discovery of transient, or “outbursting”, and “low-field” magnetars. 
Substantial progress has been made also on the theoretical side. Quite detailed 
models for explaining the magnetars’ persistent X-ray emission, the properties of the 
bursts, the flux evolution in transient sources have been developed and confronted with 
observations. New insight on neutron star asteroseismology has been gained through 
improved models of magnetar oscillations. The long-debated issue of magnetic field 
decay in neutron stars has been addressed, and its importance recognized in relation 
to the evolution of magnetars and to the links among magnetars and other families of 
isolated neutron stars. The aim of this paper is to present a comprehensive overview 
in which the observational results are discussed in the light of the most up-to-date 
theoretical models and their implications. This addresses not only the particular case 
of magnetar sources, but the more fundamental issue of how physics in strong magnetic 
fields can be constrained by the observations of these unique sources. 
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Neutron stars (NSs), the endpoint of the evolution of massive stars with 10 < M/Mq < 
25, are extremely compact remnants endowed with strong magnetic helds. Isolated (i.e. 
not belonging to a binary system) neutron stars were for a long time identihed with radio¬ 
pulsars, and only in the last two decades, mainly thanks to high-energy observations, 
was the existence of other manifestations of isolated neutron stars recognized. Among 
them, there are two groups of X-ray pulsars with remarkably peculiar properties, the 
Soft Gamma Repeaters (SGRs) and the Anomalous X-ray Pulsars (AXPs, see e.g. 
Mereghetti, 2008, for a review). The separation into two classes reflects the way in which 
these sources were originally discovered. SGRs were revealed through the detection of 
short, intense bursts in the hard X-/soft gamma-ray range (Mazets et ah, 1979b,a), and 
because of this were initially associated with gamma ray bursts (GRBs); however, burst 
emission from SGRs was soon recognized to repeat, at variance to what was observed in 
GRBs, setting the two phenomena apart. On the other hand, AXPs were identihed as 
X-ray pulsar in the soft X-ray range (< 10 keV, Mereghetti & Stella, 1995). They were 
dubbed “anomalous” because their high X-ray luminosity (~ 10^^ — 10^® erg/s) cannot 
be easily explained in terms of the conventional processes which apply to other classes of 
pulsars, i.e. accretion from a binary companion or injection of rotational energy in the 
pulsar wind/magnetosphere. Over the last few years, observations have revealed many 
similarities between these two classes of objects (see e.g. Woods & Thompson, 2006), 
including the discovery that AXPs too emit short, SGR-like bursts (Kaspi, 2000, 2003) 
and nowadays the idea that SGRs and AXPs belong to a single, unihed class is widely 
accepted. 

The main observational characteristics of SGRs and AXPs are: a) lack of evidence 
of binary companions; b) persistent (i.e. non-bursting), often variable X-ray luminosity 
in the range ~ lO^^-lO^® erg/s, emitted in the soft (0.5-10 keV) and hard (20-100 keV) 
X-ray range; c) pulsations at relatively long spin periods, clustered in the range ~ 2-12 s; 
d) large secular spin-down rate, P ~ 10“^^ — 10“^^ s/s, which, if interpreted in terms 
of electromagnetic losses from a rotating dipole in vacuo, leads to huge magnetic helds, 
~ lO^^-lO^^ G. SGRs (and AXPs, to a somewhat lesser extent) exhibit spectacular and 
frequent bursting activity, which is observed in the X-/gamma-rays on several timescales, 
ranging from snb-s to several tens of seconds. In particular, three different kinds of 
bursting events have been observed (see Sec. 5): 

• short bursts: these are the most common, with typical dnration of ~ 0.1-1 s, peak 
Inminosity of ~ lO^^-lO*^^ erg/s, and soft (~ 10 keV), thermal spectra; they are 
detected from both SGRs and AXPs; 

• intermediate bursts, which last ~ 1-40 s and have a peak Inminosity of ~ 
10"^^ — 10*^^ erg/s. These are characterized by an abrupt onset and usually also 
show thermal spectra; again, they were seen in both SGRs and AXPs; 

• giant hares. These are exceptional, rare events, with an energy ontpnt of ~ 
10"^4 _ ]^g47 erg/s, only exceeded by blazars and GRBs. They have been observed 
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only in SGRs and only three times since SGRs were discovered: from SGR 0526-66 
in 1979 (Mazets et ah, 1979a), from SGR 1900-1-14 in 1998 (Hurley et ah, 1999), 
and from SGR 1806-20 in 2004 (e.g. Hurley et ah, 2005; Palmer et ah, 2005). All 
three events started with an initial spike of ~ 0.1-0.2 s duration, followed by a long 
pulsating tail (lasting a few hundred seconds) modulated at the neutron star spin 
period. 

All together, these properties hud an explanation in the so-called “magnetar” 
scenario (Duncan & Thompson, 1992; Thompson & Duncan, 1993, 1995), according 
to which the relatively high X-ray luminosity and the bursting/outbursting activity are 
powered by the dissipation and decay of a superstrong magnetic field, ~ 10^"^-10^® G on 
the surface and possibly higher in the star’s interior. Despite no indisputable measure 
of an ultra-high magnetic held has been obtained as yet, a number of indipendent 
arguments strongly support the idea that SGRs/AXPs are indeed magnetically-powered, 
as hrst discussed by Thompson & Duncan (1995). A few of them are 

• the rotational energy loss rate E (which is believed to fuel standard pulsars) is well 
below the persistent X-ray luminosity, E Lx] 

• long spin periods (~ 10 s) can be attained in 10^ — 10^ yrs (the source age as 
inferred from that of the associated SNR) via magneto-dipolar braking only for 
helds > 10^^ G; 

• huge spin-down rates have been indeed measured in these sources, implying dipole 
magnetic helds in the range ~ 10^^ — 10^^ G; 

• the decrease of the scattering opacity in a superstrong magnetic held {B > 10^"^ 
G) pushes upwards the Eddington limit and allows a much larger luminosity to 
escape from a (magnetically) conhned plasma: this can explain the apparently 
super-Eddington luminosity of a number of bursts; 

• no stellar companions have been discovered in SGRs/AXPs, ruling out accretion as 
a possible source of energy; 

• if no more than a fraction of the magnetic energy was released in a giant hare, this 
requires B > 10^^ G; in order to power ^ 100 giant hares like that emitted in 2004 
by SGR 1806-20 over the source lifetime an internal held ~ 10^® G is needed (Stella 
et ah, 2005). 

Although alternative interpretations have been proposed (see e.g. Turolla & Esposito, 
2013, for a review and references therein), the magnetar model more naturally explains 
the properties of SGRs and AXPs, including the bursting activity and the hyper- 
energetic giant hares, and will be the focus of this review. 

Even the “persistent” emission of these sources is far from being steady. Magnetars’ 
spin-down is quite irregular, and often accompanied by glitches and timing noise. Long 
term variations in magnetars’ emission can occur either as gradual and moderate changes 
in the hux, accompanied by variations in the spectrum, pulse prohle, and spin-down rate, 
or as sudden outbursts, i.e. events during which the hux raises up to a factor ~ 1000 



Magnetars: the physics behind observations 


4 


and then decays back to a level compatible with the quiescent state over a time scale 
of months/years (see Sec. 4). Within the magnetar scenario, the hrst kind of variability 
is thought to be driven by plastic deformations in the crust which, in turn, induce 
changes in the magnetic current conhgurations. The more violent outbursts, as well as 
the glitches, the bursting activity and even the hyper-energetic giant flares could instead 
be due to sudden reconhgurations of the magnetosphere, when unstable conditions are 
reached. This may lead to crustal fractures (starquakes) and/or instabilities in the outer 
magnetosphere (possibly involving magnetic reconnection). 

Although originally discovered in the X-/soft gamma-rays, magnetars have 
been detected at different wavelengths, revealing a rich phenomenology across the 
electromagnetic spectrum. AXPs and SGRs have been discovered to emit in the optical 
and/or near-infrared (NIR) bands (e.g. Hulleman et ah, 2000; Israel et ah, 2004; Durant 
& van Kerkwijk, 2006). The optical/NIR counterparts are faint {K ~ 20) and the flux 
is only a small fraction of the bolometric flux, but still its detection can place important 
constraints on models. Several AXPs have exhibited long-term variability both in their 
optical/infrared emission and in X-rays (Israel et ah, 2002; Hulleman et ah, 2004; Rea 
et ah, 2005). Unavoidably this introduces additional uncertainties in the modelling of 
broad band spectra, based on observations at different wavelengths taken at different 
times. 

Magnetars were traditionally considered to be radio-silent, but this picture was 
challenged by the (unexpected) discovery of a pulsed radio counterpart in some sources, 
a property that seems to be peculiar to transient magnetars (Gelfand & Gaensler, 2007, 
see also Sec. 4.3 and references therein). When detected, the radio emission of magnetars 
appears to be different from that of standard radio-pulsars: the spectrum is flatter 
and the flux and pulse prohle show strong variations with time, indicating that the 
mechanisms causing the emission (or the topology of the emission region) may differ in 
the two kinds of sources. 

Association with supernova remnants or, possibly, young stellar clusters has been 
proposed for a number of sources (see e.g. Table 1 in Mereghetti, 2008; Muno et ah, 
2006; Vrba et ah, 2000; Eikenberry et ah, 2001; Figer et ah, 2005; Klose et ah, 2004), 
which, if conhrmed, leads in some cases to a progenitor with high mass (> 20 Mq), high 
metallicity, and to a relatively young age ~ 10^ yr for the neutron star (see Sec. 2.1). 

The magnetar paradigm that bursting activity is necessarily associated with a 
high dipolar held has been revolutionized by the recent discovery of a few full-hedged 
magnetars (i.e. neutron stars that displayed bursting, SGR-like activity) with a dipolar 
magnetic held comparable with that of standard radio pulsars (see Sec. 4 and references 
therein). The properties of these sources are compatible with those expected from aged 
magnetars, which may still retain a large toroidal held in the interior, occasionally 
capable of cracking the star’s crust. This discovery suggests that magnetars could be 
far more numerous than previously expected (Rea et ah, 2010; Tiengo et ah, 2013), and 
has had a number of profound implications, e.g. for star formation, supernovae, gamma 
ray bursts (see Rea, 2014b, for a complete discussion). 
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Despite the wide interest in the astrophysical community, review papers on 
magnetars were comparatively few, and mainly devoted to the diverse aspects of 
their phenomenology. Theoretical results are often scattered across many specialized 
papers, the comparison and interpretation of which are quite a challenge even for an 
informed reader. It is outside the scope of this paper to provide a detailed summary 
of magnetars’ observational properties, about which excellent review papers have been 
already published (Woods & Thompson, 2006; Kaspi, 2007; Mereghetti, 2008; Hurley, 
2011b; Rea & Esposito, 2011); an updated list of sources, containing all the essential 
data, is available in the online McGill magnetar cataloguel (Olausen & Kaspi, 2014), 
and while preparing this review we also created a Magnetar Burst Library which is 
now maintained by the Univ. of Amsterdam §. Here we will focus on the theoretical 
interpretation of the emission properties of magnetars and on a cross comparison of 
the models presented so far to describe them. A brief summary of the observational 
properties, which is not necessarily complete but sets the context for the subsequent 
discussion, is placed at the beginning of each section, when needed. Our main aims are 
to review the state of the art in the theoretical modelling, to outline which observational 
facts are robustly explained by current models and to discuss the open issues which still 
remain to be addressed. 

The paper is organized as follows. We begin with a summary of the mechanisms that 
can lead to the birth of a highly magnetized neutron star, discuss how magnetars evolve 
and briefly touch the link between magnetars and other classes of Galactic, isolated 
neutron stars (Sec. 2). Sec. 3 is dedicated to the twisted magnetosphere model and its 
ability to explain the observed persistent emission in different wavebands. Transient 
magnetars and their observations in the radio band are reviewed in Sec. 4, while Sec. 5 
contains a thorough discussion of burst emission and magnetar seismology. Gonclusions 
follow. 

2. Birth and evolution of a magnetar 

2.1. Magnetars formation 

According to the original picture by Duncan and Thompson (Duncan & Thompson, 
1992; Thompson & Duncan, 1993), ultra-magnetized neutron stars form through 
magnetic held amplihcation by a vigorous dynamo action in the early, highly convective 
stages. Rotation and convection produce two types of dynamo effects in an astrophysical 
plasma: the a dynamo, arising from the coupling of convective motions and rotation, 
and the u dynamo, driven by differential rotation. In proto NSs both effects are present 
and since the a-u dynamo operates at low Rossby numbers, the initial spin period 
must be short, < 3 ms, to ensure efficient convective mixing (Duncan & Thompson, 
1992). Magnetars would be, then, the endpoint of the evolution of massive stars 

f The on line McGill catalogue can be found at 

http: //WWW.physics.mcgill.ca/ ~ pulsar/magnetar/main.html. 

§ See the Amsterdam Magnetar Burst Library, http://staff.fnwi.uva.n1/a.l.watts/magnetar/mb.html 
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with rapidly rotating cores. Rapidly spinning, collapsing stellar cores are expected 
to produce highly energetic supernovae (Duncan & Thompson, 1992; Thompson et ah, 
2004; Bucciantini et ah, 2007), because a signihcant fraction of the rotational energy, 
Erot ~ 3 X 10^^(P/1 ms)“^ erg, is transferred to the ejecta via the strong magnetic 
coupling with the proto-neutron star. Any observational signatures that magnetars were 
born in (above-average) energetic events were searched for in a number of supernova 
remnants positively associated with SGRs/AXPs, but no evidence has yet been found 
(Vink & Kuiper, 2006; Vink, 2008). If the internal magnetic held is ~ 10^® G, however, 
rotational energy can be efficiently carried away by gravitational waves, which do not 
interact with the ejecta (Dall’Osso, Shore & Stella, 2009). 

Alternatively, it has been suggested that ultra-strong helds in neutron stars result 
from magnetic hux conservation (the fossil held scenario; Ferrario & Wickramasinghe, 
2006, 2008). Ferrario & Wickramasinghe (2006), starting from a parameterized model 
of the distribution of magnetic hux on the main sequence and of the spin period of 
neutron stars at birth, derived the expected properties of isolated radio pulsars in the 
Galaxy, given the spatial distribution of the initial mass function and star formation 
rate. Gomparison with the data in the 1374-MHz Parkes Multi-Beam Survey was then 
used to constrain the model parameters. They hud that the distribution of the magnetic 
held in the core of the OB progenitors comprises ~ 8% of stars with a magnetic held in 
excess of ~ 1000 G. The core-collapse supernovae of these high-held stars can produce 
~ 25 magnetars with properties (surface magnetic held, spin period, age) in agreement 
with those observed in SGRs/AXPs. As hrst noted by Spruit (2008), the number of 
Galactic magnetars predicted by the fossil held model may be too low, and this is a 
more and more serious issue, given the steady increase of the magnetar population and 
the possibility that many “dormant” SGRs/AXPs lurk among “standard” radio pulsar 
(see Sec. 4). A possibility is that magnetars are formed through diherent channels: 
for instance it has been suggested that at least part of the magnetars may be born as 
rapidly rotating neutron stars in systems in which the core of the collapsing star was 
accelerated by tidal synchronization in a very close binary (Popov & Prokhorov, 2006; 
Bogomazov & Popov, 2009). 

Interestingly, the high-held progenitors of magnetars should be in the far end of the 
mass distribution of OB stars, with masses ~ 20-45 Mq, which, in standard evolutionary 
models, should mostly have given rise to black holes (Ferrario & Wickramasinghe, 2008, 
see also Glark et al. 2005). The notion that magnetars descend from massive stars 
(typically above the canonical neutron star-black hole divide) received further support 
from the observational evidence that (some) SGRs/AXPs are associated with young 
clusters of massive stars. The progenitor mass of SGR 1806-20 and the AXP IE 1048.1- 
5937 has been estimated to be in excess of ~ 30 Mq (Bibby et ah, 2008; Gaensler et 
ah, 2005a); the progenitor mass of SGR 1900-1-14 appears, however, to be ~ 17 Mq 
(G lark et ah, 2008; Davies et ah, 2009). One of the strongest evidence in favour of 
high-mass progenitors of SGRs/AXPs came from the robust association of the AXP 
GXO J164710.2-455216 with the young cluster Westerlund 1 (Muno et ah, 2006). Since 
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the cluster is ~ 4 Myr old (Clark et ah, 2005), the minimum mass of a star that could 
have reached the supernova stage is ~ 4OM0. Hence the claim that CXO J164710.2- 
455216 originated from a star with M > 40Mq. Very recently Clark et al. (2014) 
proposed that CXO J164710.2-455216 was born in a massive binary and found evidence 
for the former companion, the runaway star Wdl-5, ejected from the system when 
the magnetar progenitor exploded. If Wdl-5 and CXO J164710.2-455216 were indeed 
related, evolution in the binary would lead to a decrease of the progenitor mass through 
strong mass loss when it entered a Wolf-Rayet phase, and common envelope evolution 
would prevent spin-down of its core. Magnetar birth in a binary may then be a key 
ingredient to bring the progenitor mass within the neutron star formation range, and to 
provide the high core rotational speed required for the onset of the convective dynamo. 

Magnetars are also increasingly popular as the central engine powering GRBs, 
following the original suggestion by Usov (1992, see also Zhang & Meszaros 2001; 
Metzger et al. 2011). Ultra-magnetized neutron stars have been invoked to explain 
the properties of both short and long GRBs. The proto-magnetar would result from 
coalescence in a double-degenerate binary (or accretion-induced collapse of a white 
dwarf) in this hrst case and in a core-collapse supernova in the second (e.g. Paczyinski, 
1986; Rosswog, Ramirez-Ruiz & Davies, 2003; Giacomazzo & Perna, 2013; Metzger, 
Quataert & Thompson, 2008; Woosley, 1993; MacFayden & Woosley, 1999). Indeed, a 
signihcant fraction of the Swift long GRBs exhibit late flares and plateau phases in the 
lightcurve that provide evidence for longevity and on-going activity of the central engine 
(see e.g. Curran et al., 2008; Margutti et al., 2010; Bernardini, et al., 2011b; Nousek 
et ah, 2006; Zhang et al., 2006). The plateau, which occurs around 10^ — 10"^ s after 
the trigger, has a fluence that can be as high as the fluence of the prompt emission. 
According to the magnetar model the plateau phase is powered by the initial spin-down 
of a newly born magnetar, powering a relativistic wind (Fan et al., 2006; Troja et al., 
2007; Lyons et al., 2010; Dall’Osso et al., 2011; Bernardini et ah, 2012; Metzger et al., 
2011). Moreover, Rowlinson et al. (2013) have recently shown that 18 of the Swift short 
GRBs (i.e. 64% of the entire sample) can be clearly htted with a magnetar plateau 
phase, while for the rest the quality of the data is insufficient to prove or exclude the 
presence of the plateau. Out of 18 robust candidates, 10 are thought to collapse later 
to a black hole, while the others may have left behind a rapidly rotating new magnetar. 
Although these studies are not a direct, conclusive proof of the magnetar paradigm, 
they certainly indicate the frequent occurrence of late central activity, which has crucial 
implications for the origin of the central engine. A smoking gun that may allow to 
differentiate between models would be the detection of gravitational waves associated 
to the event (Rowlinson et ah, 2013, and references therein). 

Another link between magnetars and GRBs has been proposed following the 
observations of giant flares. Since all these events started with an initial, very energetic 
sub-s spike, it has been proposed that giant flares, if emitted by extragalactic SGRs, 
may appear at Earth as short gamma-ray bursts (Palmer et ah, 2005; Hurley et al., 2005; 
Hurley, 2011a). The main causes of uncertainty for proving this idea are in maximum 
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energy released in the flare and in the spectral properties of the narrow peak. By 
considering the flare emitted by SGR 1806-20 and by varying the assumptions about 
the peak spectral shape, Popov & Stern (2006) computed the possibility of detection by 
BATSE of giant flares with an energy of 10^^ or 10^® erg, as a function of the distance. 
They found that the first kind of event can be seen up to a few Mpc (therefore in M82, 
M83, NGC253 and NGG4945), while the second can in principle be visible up to the 
Virgo cluster. However, as already noted by Popov & Stern (2006), this prediction may 
be too optimistic, since no evidence has been found for an excess of BATSE short GRBs 
from the direction of M82, M83, NGG253 and NGG4945, nor from the Virgo cluster 
(Palmer et ah, 2005). Similarly, negative results have been reported by Lazzati et ah 
(2005); Tanvir et al. (2005), and overall these studies suggest that no more than a few 
percent, maybe up to ~ 8% of the short GRBs seen by BATSE could be giant flares 
from extragalactic SGRs (see also Hurley et ah, 2005; Growther et ah, 2011; Svinkin et 
ah, 2015, the latter for a recent update on the detection upper limits). 

2.2. Magneto-thermal evolution 

A major issue in establishing the magnetic evolution of NSs (and of magnetars in 
particular) is that observations place very little, if any, constraint on the structure and 
strength of the internal magnetic held. Glearly, in a magnetar the internal held must 
be strong enough to sustain the source activity and its geometry must allow magnetic 
energy to be released. While there are several indications that the large-scale, external 
held can be reasonably assumed to be (nearly) dipolar, the internal held most likely 
contains both toroidal and poloidal components (e.g. Geppert, Kiiker & Page, 2004, 
2006, and Sec. 3.1). A further complication comes from the at present poor knowledge 
of where the internal held resides. The held can either permeate the entire star (“core” 
fields), or be mostly confined in the crust (“crustal” fields), depending on where its 
supporting (super)currents are located. 

The more general configuration for the internal held in a NS will be, then, that 
produced by the superposition of current systems in the core and the crust. As stressed 
by Pons & Geppert (2007), the relative contribution of the core/crustal fields is likely 
different in different types of NSs. In old isolated radio pulsars, where no held decay is 
observed, the long-lived core component may dominate, while a sizeable, more volatile 
crustal held is probably present in magnetars, for which substantial held decay over a 
timescale ~ 10^-10^ yr is expected (e.g. Goldreich & Reisenegger, 1992). As pointed 
out by Glampedakis, Jones & Samuelsson (2011) ambipolar diffusion plays little role 
in magnetar cores during their active lifetime (after crystallization, the absence of 
convective motions already quenches ambipolar diffusion in the crust). Therefore, if 
the decay/evolution of the magnetic held is indeed the cause of magnetar activity, it 
is likely to take place outside the core and be governed by Hall/Ohmic diffusion in the 
stellar crust. Other mechanisms, e.g. flux expulsion from the superconducting core, 
due to the interaction between neutron vortices and magnetic flux tubes, are highly 
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uncertain and very difficult to model. For these reasons, recent investigations of the 
magnetic held evolution in magnetars has focused only on the crustal component of the 
held. 

The relative importance of the Ohmnic decay and Hall drift is strongly density- and 
temperature-dependent. Thus, any self-consistent study of the magnetic held evolution 
must be coupled to a detailed modelling of the neutron star thermal evolution, and vice 
versa. This basically means that the induction eqnation for B must be solved together 
with the cooling, a qnite challenging nnmerical task. Early ehorts in this direction nsed 
a split approach. Pons & Geppert (2007) stndied the evolntion of the held by solving 
the complete indnction eqnation in an isothermal crnst, bnt assnming a prescribed time 
dependence for the temperature. They found that crustal magnetic helds in NSs snher 
signihcant decay dnring the hrst ~ 10® yr and that the Hall drift, althongh inherently 
conservative (i.e. alone it cannot dissipate magnetic energy), plays an important role 
since it may reorganize the held from the larger to the smaller (spatial) scales where 
Ohmic dissipation proceeds faster. 

The cooling of magnetized NSs with held decay was investigated by Agnilera et ah 
(2008, see also Agnilera et ah 2009; Kaminker et ah 2006, 2007, 2009) by adopting 
a simple, analytical law for the time variation of the held which incorporates the 
main featnres of the Ohmic and Hall processes. The fnlly conpled magneto-thermal 
evolution of a NS was addressed by Pons et ah (2009), inclnding all realistic microphysics. 
However, owing to nnmerical difhcnlties in treating the Hall term, their models acconnt 
only for Ohmic dihusion. A complete treatment of magneto-thermal evolntion, properly 
inclnding the Hall term, was recently presented by Vigano et al. (2013, see also Vigano, 
Pons & Miralles 2011b). Their calculations conhrm the basic pictnre ontlined in Pons et 
al. (2009), althongh the presence of the Hall drift introdnces some remarkable diherences. 
Contrary to the pnrely dissipative case, evolution is not very sensitive to the initial 
relative strength of the toroidal component with respect to the poloidal one, nnless the 
former is mnch higher than the latter. This is because a toroidal component bnilds np 
anyway dne to the Hall term, even starting with a pnrely poloidal held conhgnration. 
Models are not strongly dependent on other parameters (notably the star mass) either, 
so that the evolntion is mostly controlled by the initial valne of the dipolar held. Fig. 1 
shows the evolntion of the dipolar held and of the thermal Inminosity for diherent initial 
magnetic geometries: core held (model B14), core-|-crustal held (C14) and pnrely crustal 
held (A14, AT14). The diherent decay pattern of crustal vs. core helds is evident. 

2.3. Magnetars and other neutron star classes 

Over the last two decades onr pictnre of the Galactic nentron star population has 
changed drastically, mainly thanks to high-energy observations. Besides SGRs/AXPs, 
the existence of several new classes of isolated nentron stars (INSs), with properties 
qnite at variance with those of ordinary radio-pulsars (PSRs), has emerged: the central 
compact objects in snpernova remnants (GGOs in SNRs), the rotating radio transients 
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Figure 1 Left panel: evolution of the dipole field polar strength. Right panel: same 
for the thermal luminosity. In all models the initial poloidal held is 10^^ G; the initial 
toroidal held is zero apart from model A14T, where it is 5 x 10^^ G. The star mass is 
1.4Mq (from Vigand et al., 2013, with OUP permission). 


(RRaTs) and the X-ray dim INSs (XDINSs or M7) (see e.g. Kaspi, 2010; Harding, 2013; 
De Luca, 2008; Ho, 2012; Burke-Spolaor, 2012; Turolla, 2009, for reviews). All these 
sources are radio-silent or, in the case of RRaTs (and SGRs/AXPs), show only sporadic 
(transient) radio emission (see Sec. 4.3). They were discovered as X-ray pulsators, with 
the exception of the RRaTs (only one is currently known as an X-ray source, McLaughlin 
et ah, 2007), and their spectrum is mostly thermal. While the periods are quite long 
(from ^ 0.1-0.4 s for the GGOs to ~ 1-10 s for the XDINSs and RRaTs), their period 
derivatives span a large interval, with implied magnetic helds ranging from as low as 
~ 3 X 10^° G in some of the GGOs (which are sometimes referred to as the “anti- 
magnetars”), to ~ 10^^ —10^^ G in RRaTs and XDINSs (see Keane et al., 2011; Turolla, 
2009). Ages are also very diherent, GGOs being quite young (the associated SNR age 
is < lO'^ yr) and XDINSs much older (the estimated dynamical ages are ~ 10^ yr, e.g. 
Mignani et al., 2013, and references therein); in both cases the “true” ages turn out to be 
shorter than the spin-down ages. Like PSRs, RRaTs appear to be rotationally-powered, 
while the (thermal) X-ray emission from XDINSs and GGOs is powered by the release 
of residual heat. The position of the various sources in the P P diagram is shown in 
Fig. 2. 

Although the number of detected sources in each class is fairly limited in comparison 
to that of PSRs (7 XDINSs, > 70 RRaTs, 8 GGOs, about 20 SGRs/AXPs, and few 
candidates in each class, vs. > 2000 PSRs||), the estimated birthrate of XIDNSs and 
RRaTs is comparable to and possibly higher than that of PSRs, /dpsR ~ 0.015-0.03 yr“^ 
(Popov, Turolla, & Possenti, 2006; Keane & Kramer, 2008, and references therein). The 
magnetar birthrate is lower than those of other classes, /Smag ~ 0.003 yr“^, although this 


II ATNF pulsar catalogue, http://www.atnf.csiro.au/research/pulsar/psrcat/ 
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Figure 2 The P P diagram illustrating the placement of the different isolated neutron 
star classes. The blue dots mark pulsars detected both in the radio and X-ray bands, the 
red ones those observed only at X-ray energies. The lines of constant age and magnetic 
held are also shown (courtesy R.P. Mignani). 


is likely a lower limit given the increasing number of SGRs/AXPs discovered recently. 

This clearly is an issue, since the sum of the birthrates of the various INS types 
cannot exceed the core-collapse supernova rate in the Galaxy, /3 sn ~ 0.02 ± 0.01 yr“^. 
Unless the current hgures for the INS birthrates are grossly overestimated (and/or INSs 
can form through other channels), this implies that some evolutionary links exist among 
the different classes (Keane & Kramer, 2008). That XDINSs might be aged, worn-out 
magnetars has been suggested repeatedly, on the basis of the similarity of the periods and 
the (relatively) high magnetic helds of the former (e.g. Turolla, 2009). Besides the need 
to hnd evolutionary links among the groups, the variety of INS manifestations brings 
in an even more fundamental question: which initial parameters determine whether a 
proto NS will become, say, a magnetar or a PSR ? The idea that the properties (and 
the evolution) of an INS are governed by a limited number of macrophysical quantities 
at birth (e.g. mass, magnetic held, period) may indeed open the way to what has been 
called the “grand unihcation of neutron stars”, or GUNS for short (e.g. Kaspi, 2010; 
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Igoshev, Popov & Turolla, 2014). 

Magnetic field decay is bound to play a central role in any attempt to build a GUNS. 
Popov et al. (2010) were the first to perform INS population synthesis calculations 
including magneto-thermal evolution, adopting the treatment of Pons et al. (2009). 
Their model satisfactorily reproduces all INS populations if the initial magnetic field 
follows a log-normal distribution with a mean value Bq = 1.8 x 10^^ G. Their picture 
conhrms that the magnetic held decays substantially (by a factor > 10) in the most 
magnetic stars, but provides no clear indications for evolutionary links among the 
different INS groups. New population synthesis calculations, including more updated 
magneto-thermal evolutionary models, have been recently presented by Gullon et al. 
(2014). A more decisive indication that such links indeed exist comes from the tracks 
computed by Viganb et al. (2013) by coupling the magnetic held and period evolution 
(see Fig. 3). The main ehect of magnetic held decay is to produce a sharp bending of 
the track downwards after a time ~ 10^ yr for strong initial helds. This implies that 
the star’s period does not increase indehnitely but freezes at an asymptotic value which 
depends on the initial magnetic held, the mass of the star and the crust resistivity (see 
also Dall’Osso, Granot & Piran, 2012). A comparison between the theoretical tracks 
and the positions in the P P plane of INSs of diherent types (see again Fig. 3) suggests 
that “moderate” magnetars {Bq = a few x 10^^ G) evolve into XDINSs. 

3. Persistent emission 

3.1. Magneto spheric twist 

The current picture of a magnetar magnetosphere relies on the notion that the star’s 
external magnetic held dihers from a simple, potential dipole, which is usually assumed 
to be the case for “standard” neutron stars. The reason for which the external held is 
not dipolar is to be sought in the structure of the internal magnetic held. Over the last 
decade, analytical and numerical investigations have shown that any stable conhguration 
for the internal magnetic held of a star has necessarily to contain both a poloidal and 
a toroidal component (e.g. Tayler, 1973; Flowers & Ruderman, 1977; Braithwaite & 
Spruit, 2006; Braithwaite & Nordlund, 2006; Braithwaite, 2008, 2009). In particular, 
Braithwaite (2009) investigated stable, axisymmetric magnetic equilibria and found that 
the ratio of the two components must be such that aE/U < Ep/E < 0.8, where E and 
U are the total magnetic and gravitational energies, Ep is the energy associated with the 
poloidal component and a is a numerical factor. Given that E/U < 10“^^ and a ~ 10^ for 
a neutron star, its internal magnetic held likely comprises a toroidal component at least 
of the same order as, and possibly much stronger than, the poloidal one. The instability 
of purely poloidal or toroidal magnetic conhgurations was proven also by Newtonian 
(e.g. Lander & Jones, 2011a,b) and general-relativistic (e.g. Giolh et ah, 2011; Giolh 
& Rezzolla, 2012) numerical simulations (see also Giolh, 2014, for a recent overview). 
Although, earlier attempts with the twisted-torus model (a likely conhguration for the 
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Figure 3 Evolutionary tracks in the P P plane of INSs with different initial magnetic 
fields. Asterisks mark the real age of the source along the track (10^, 10^, 10®, 5 x 10® yr) 
and the dashed lines give the tracks with constant B. MAG = SGRs/AXPs, XIN = 
XDINSs, HB = high-B PSRs, RPP = PSRs (from Vigand et ah, 2013, with OUP 
permission). 


internal stellar held, e.g. Braithwaite & Nordlund, 2006) pointed towards poloidal- 
dominated geometries, which are themselves unstable (Giolh, 2014, and references 
therein), more recent calculations indicate that large toroidal helds (comprising up 
to 90% of the total magnetic energy) can indeed be achieved in this framework 
(Giolh & Rezzolla, 2013, see also TAkgiin et ah, 2013 for magnetic conhgurations with 
Btor 3> Bpoi). Due to the complexity of the problem, most of those studies considered 
either the internal held structure (given an assumption for the magnetosphere) or the 
external magnetosphere (assuming an internal current distribution). The hrst global 
models, recently presented by Ruiz et ah (2014); Glampedakis et ah (2014), and Pili 
et ah (2015) in both Newtonian gravity and GR, appear promising, although a proper 
analysis of their stability has not been carried out yet. 

In a magnetar, where the internal held can exceed 10^® G, magnetic stresses can 
overcome the crustal tensile strength (Thompson & Duncan, 1995). The easiest way in 
which the crust reacts to the applied forces is through horizontal displacements, parallel 
to the magnetic equipotential surfaces, i.e. a magnetically-stressed crustal patch tends to 
rotate by an angle A(j) (Thompson et ah, 2000). This can be understood by considering a 
hux tube in which the toroidal component is non-zero in the crust and vanishes outside 
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Figure 4 A schematic view of a magnetar internal field (Thompson & Duncan, 2001, 
©AAS. Reproduced with permission. A link to the original article via DOI is available 
in the electronic version). 


the star (Thompson, Lyutikov & Kulkarni, 2002). Because the conductivity is much 
higher in the star’s interior, the currents supporting the non-potential R-held will close 
in a thin surface layer. The Lorentz force acting on the current, and hence on the layer, 
is Fl = jxB/c = jx {Bp -|- Bt)/c, where j is the current density. The part of 
due to the toroidal held Bt tries to produce a vertical displacement, which is unlikely 
to occur due to the strong stratihcation (Reisenegger & Goldreich, 1992), while the 
part associated with the poloidal component Bp results in a slippage in the horizontal 
direction (see Fig. 4). 

A direct effect of the magnetically-induced rotation of a surface platelet is the 
twisting of the external held. Since the external magnetic held lines are anchored to 
the crust, a torsional displacement of the surface layers produces a transfer of magnetic 
helicity from the interior to the exterior. If the external held is initially dipolar, it 
will acquire a non-zero toroidal component, a twist, conhned to the held lines whose 
footpoints are on the displaced layer. In a twisted magnetosphere, currents necessarily 
how also along the closed held lines to support the non-potential held. This is at 
variance with what is usually assumed to occur in “normal” radio-pulsars, where charges 
(the Goldreich-Julian currents) move only along the open held lines (again because 
the R-held is non-potential in that region). The presence of large-scale currents in 
a magnetar magnetosphere has major implications in shaping the emergent spectrum 
through repeated resonant cyclotron scatterings, as will be discussed in Sec. 3.3. The 
gradual implant and subsequent decay of a magnetospheric twist has been often invoked 
to explain the long terms evolutions of some magnetars (Mereghetti et ah, 2005b; 
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Campana et al., 2007), for example the behaviour observed before or after a series 
of bursts or a giant flare. For instance, before the giant flare emitted by SGR 1806-20, 
the source properties changed remarkably: a study of the observed long-term variations 
indicated a clear correlation among the increases in spectral hardening, spin-down rate, 
and bursting activity (Thompson, Lyutikov & Kulkarni, 2002; Mereghetti et ah, 2005b). 
The proposed scenario assumes the onset of a gradually increasing twist: this, in fact, 
results in an increasing optical depth for resonant cyclotron scattering, and causes a 
progressive hardening of the X-ray spectrum. At the same time, the spin-down rate 
is expected to increase because, for a fixed dipole field, the fraction of held lines that 
open out across the speed-of-light cylinder grows. Since both the spectral hardening 
and the spin-down rate increase with the twist, the model predicts that they should be 
correlated in agreement with the observations. 

Although magnetospheric twists are expected to be localized, meaning that they 
do not affect the entire magnetosphere (Thompson, Lyutikov & Kulkarni, 2002; 
Beloborodov, 2009), nearly all studies on the properties of the persistent emission from 
magnetars rely on the “globally twisted magnetosphere” hrst proposed by Thompson, 
Lyutikov & Kulkarni (2002). In this model it is assumed that the external magnetic 
held is initially dipolar and that, as a consequence of crustal displacements, a certain 
amount of shear is added to the held. If one restricts to magnetostatic equilibria in a 
low-density plasma, the momentum equation reduces to^ j x B = 0, which, combined 
with the Ampere-Maxwell equation V x B = (47r/c)j gives the force-free condition 

{V X B) X B = 0. (1) 


By expressing the poloidal component through the hux function V, an axisymmetric 
held has the most general form 


B 


VV{r,9) X 
r sin 6 


+ B^{r, 9)u^ , 


( 2 ) 


where B^ is the toroidal component and the unit vector in the 0 direction. By 
exploiting the force-free condition one can explicitly write the magnetic held as 


B 



f, vl_ [cV 

’ sin 6* y p -|- 1 sin 9 


( 3 ) 


where a prime denotes a derivative with respect to /i = cos 9, Bp is the polar value of 
the magnetic held, R^s is the star radius, G is a constant and 0 < p < 1 is the radial 
index. The function /(p) satishes the Grad-Shafranov equation 


(l-p2)/" + p(p + i)/ + c/i+2/P = o 


( 4 ) 


which is a second order ordinary diherential equation for the angular part of the hux 
function. Since equation (4) must be (numerically) solved subject to three boundary 
conditions (Thompson, Lyutikov & Kulkarni, 2002; Pavan et ah, 2009), the constant 

^ SGRs/AXPs are slow rotators and the Coulomb force is negligible in the inner magnetosphere. 
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Figure 5 A globally twisted dipolar field (right panel) as compared with a pure dipole 
magnetic configuration (left panel). 

C is an eigenvalue and is completely specified once p is fixed. The solution of 
equation (4) completely determines the external magnetic field and provides a sequence 
of magnetostatic, globally-twisted dipole fields by varying the index p. A picture 
illustrating a globally-twisted dipole is shown in Fig. 5. 

Besides controlling the radial decay, the value of p also fixes the amount of shear of 
the field. In fact, the twist angle, i.e. the angle through which a field line has rotated 
when it comes back to the stellar surface, is defined as 




The effect of decreasing p is to increase with respect to the other components, and 
consequently to increase the shear. The limiting values p = 0,1 correspond to a split 
monopole and an untwisted dipole, respectively. 

Primarily to assess the role played by the magnetic geometry on the emergent 
spectra, the effects on the spectra of other sheared magnetospheric configurations have 
been investigated. In these models the helicity is not uniformly distributed and, in a 
sense, they can be thought of as closer to the realistic case in which the twist is localized. 
Globally-twisted multipolar fields have been considered by Pavan et ah (2009), following 
essentially the same approach adopted by Thompson, Lyutikov & Kulkarni (2002) for 
the dipole. More recently, Vigano et al. (2012, see also Viganb, Pons & Miralles 2011) 
explored more general, non self-similar, force-free configurations for the external 5-field 
in which an arbitrary function is used to control the spatial distribution of the twist. 
The implications for spectral calculations will be discussed in Sec. 3.3.2 

Once implanted by a starquake, the twist must necessarily to decay. In a genuinely 
static twist {dB/dt = 0), in fact, the electric and magnetic fields are orthogonal. This 
implies that the voltage drop between the footpoints of a field line vanishes since = 0, 
so that there is no force that can extract particles from the surface and lift them against 
gravity thus initiating the current required to sustain the sheared field, Jb = cV x B/Att. 
As discussed by Beloborodov & Thompson (2007, see also Thompson et al. 2000), the 








Magnetars: the physics behind observations 


17 


twist decays precisely to provide the potential drop required to accelerate the charges. 
A non-vanishing is maintained by self-induction and the twist evolution is regulated 
by the balance between the conduction current j and Jb, dE\\/dt = If j < js 

the magnetosphere becomes charge starved and E\\ grows at the expense of the magnetic 
held, injecting more charges into the magnetosphere. On the other hand, when j > Jb 
the held decreases, reducing the current. The magnetosphere is then in a dynamical 
(quasi-)equilibrium with j ~ Jb over a time-scale < fdecay The potential drop across a 
held line is maintained close to the pair production threshold, e<h ~ 1 GeV, and the rate 
of magnetic energy dissipation is E^ag ~ /*h, where I ~ JbI^ is the current and I the 
linear size of the twisted region (Beloborodov & Thompson, 2007). The magnetic energy 
stored in the twist is Emag ~ I'^Rns/c^, and the twist decay time tdecay ~ Emag/Emag 
turns out to be ~ 1 yr for typical parameter values. The detailed evolution of a twisted 
magnetosphere has been investigated by Beloborodov (2009). 

3.2. Current distribution 

A twisted magnetosphere can be regarded as a force-free conhguration, threaded by 
currents that how along the 5-held lines with j ~ Jb- Charges are extracted from 
the star’s surface and accelerated by the electric held parallel to B. In the simplest 
picture (Thompson, Lyutikov & Kulkarni, 2002), the charge how consists of two counter¬ 
streaming currents: electrons and ions moving in opposite directions, so that charge 
neutrality is ensured. Using a simple, unidimensional circuit analogue, a twisted hux 
tube is akin to a relativistic double layer (Beloborodov & Thompson, 2007; Carlqvist, 
1982), in which electrons/ions leave the anode/cathode and are accelerated by a 
potential drop, which, in turn, depends on the current. Beloborodov & Thompson 
(2007) pointed out that such a conhguration cannot be realized in the magnetosphere 
of a magnetar. In order to produce j ~ j_b, in fact, the Lorentz factor of the electrons 
needs to be sufficiently high (7 ^ 10 ®) to make one-photon pair production in the 
strong magnetic held through resonant cyclotron up-scattering unavoidable well before 
7 attains such large values. Currents are expected to be carried mostly by pairs, the 
corona being in a state of self-organized criticality with a voltage drop near the threshold 
for the ignition of pair cascades. 

The analysis by Beloborodov & Thompson (2007) revealed much of the (complex) 
physics of a twisted magnetosphere. Still, being based on an idealized circuital model, 
it was not particularly suited for being used in spectral modelling. For this reason 
most investigations in this direction have resorted to the simpler, albeit less physically 
sound, picture of electron/ion currents. Under this assumption and having specihed 
the magnetic conhguration, the density of magnetospheric particles is automatically 
hxed once the particle velocity is known. In particular, for a force-free globally twisted 
dipolar held (Thompson, Lyutikov & Kulkarni, 2002; Fernandez & Thompson, 2007; 
Nobili, Turolla & Zane, 2008a), the charge density can be derived from the condition 
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j = 3 b 


ne(f,/3) 


P+1 / B 
Ane \Be) r|(/3)| ’ 


( 6 ) 


where {(A) is the average charge velocity (in units of c). The previous expression gives 
the co-rotation charge density of the space charge-limited flow of ions and electrons 
from the NS surface, that, due to the presence of closed loops in a twisted held, is 
much larger than the Goldreich-Julian density, ugj- In a general scenario, positive 
and negative charges (with densities n±) how in opposite directions with velocities v± 
and j = is = e (n+n+ — n_n_), where < 0. Electrons are assumed to how from 

north to south and conversely for ions. This breaks the symmetry between the star’s 
two hemispheres, and, for instance, implies that the observed spectrum will be diherent 
when viewed from the north or the south pole (see Nobili, Turolla & Zane, 2008a). The 
presence of ions introduces negligible ehects on the continuum spectra. Photons may 
still scatter oh ions, which are heavier and concentrated toward the star’s surface, but 
this is likely to give rise at most to a narrow absorption feature at the ion cyclotron 
energy (Thompson, Lyutikov & Kulkarni, 2002; Fernandez & Thompson, 2007; Tiengo et 
ah, 2013, and discussion in Sec. 4.2). For this reason, these models are often referred to 
as “unidirectional hows”, with reference to electrons only, while the term “bidirectional 
hows” is used when pairs are accounted for. 

In the absence of any detailed modelling of the current how (e.g. particle 
acceleration, interaction of charges with radiation traversing the magnetosphere), the 
velocity distribution is assumed to be spatially independent so that the charge velocity 
is a free parameter of the model. A major diherence between the various models 
(Fernandez & Thompson, 2007; Nobili, Turolla & Zane, 2008a) is in the adopted 
description of the velocity distribution of the scattering particles. In a strong magnetic 
held the electron distribution is expected to be largely anisotropic: e~ stream freely 
along the held lines, while they are conhned in a set of cylindrical Landau levels in the 
plane perpendicular to B. For this reason, Nobili, Turolla & Zane (2008a) assumed 
a collective (bulk) electron motion with velocity v^uik associated with the charge how 
in the magnetosphere, superimposed on a 1-D relativistic Maxwellian distribution at a 
given temperature Te which simulates the particle velocity spread (and the dissipation 
due to local turbulence and possible instabilities). 

The (invariant) distribution function is then 


dUe 

d{^/3) 


He exp (- 770 e) 
2iLi(l/0,) 




(7) 


where 7 ' = 776 «zfc(l — /d/Sbuik), ©e = kTe/nieC^, Ki is the modihed Bessel Function of the 
hrst order and fe = 'y~^n~^dne/d/3 is the momentum distribution function. 

In this model electrons are, then, assumed to move isothermally along the held 
lines, whilst at the same time receiving the same boost from the electric held. This is 
at least in qualitative agreement with the results of the simplihed bidirectional model 
by Beloborodov & Thompson (2007). A diherent choice was made by Fernandez & 
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Thompson (2007), who did not include the charge bulk motion in their models (despite 
this being a necessary ingredient to reproduce the current flow) and assessed the effects of 
other possible (local) distributions, either thermal or not thermal, in a few representative 
cases. In particular, they considered: 

a) a mildly relativistic, 1-dimensional flow described by a Boltzmann distribution at a 
temperature /crTo = (bo — l)meC^ 

1 [ ^ ' 

70-ij ’ 

b) a mildly relativistic, one-dimensional gas with the same Boltzmann distribution but 
extending over positive and negative momenta, in order to simulate an electron-positron 
flow; and 

c) a broad power-law in momentum, 

fi/3'y) oc (/^y)", (9) 



which mimics a warm relativistic plasma. 

As we discuss in the next section, charges must flow at mildly relativistic speed (7 ~ 
1 ) in the twisted magnetosphere for the model to successfully reproduce the observed 
X-ray spectra. While this is not a problem for the (over) simplified unidirectional flows 
discussed earlier on, where the velocity is a tunable parameter, the question of what 
occurs in a more realistic description which includes pairs is a crucial one. As discussed 
by Beloborodov & Thompson (2007), in a twisted magnetosphere electrons and ions, 
lifted from the star’s surface and accelerated by the self-induction electric held, must 
efficiently produce e^, at least if the current circulating in the circuit is ~ js. According 
to their analysis, flow with highly relativistic speed (7 ~ 10 ^) and a large velocity 
spread in the inner magnetosphere (r ~ Rns)- This poses a problem for the mildly- 
relativistic, counter-streaming model which is only valid in the (unphysical) assumption 
that pair production is neglected (Beloborodov, 2013a). On the other hand, in the 
presence of pairs, the electric field along the 5-hnes, E\\, is incapable of counteracting 
the radiative pull outwards because, at the same time, it acts as an accelerator for 
the charges of opposite sign. The result is that charges are accelerated outwards at 
relativistic velocity and no self-consistent solution yielding mildly relativistic flows is 
possible. 

Pair production in a twisted magnetosphere has been investigated in several works. 
As discussed by Medin & Lai (2007), for an iron crust and magnetic fields as high 
as ~ lO^^G, vacuum gaps may be formed above the polar regions of SGRs/AXPs, 
with subsequent pair creation. Near the stellar surface, where the magnetic held B 
exceeds the quantum limit Bq ~ 4.4 x 10^^ G, scattering between fast electrons and 
~ 1 keV seed photons generates high-energy gamma rays that immediately convert 
to electron/positron pairs via one-photon pair production (e.g. Harding & Lai 2006). 
This idea, originally proposed by Beloborodov & Thompson (2007), has been more 
recently reconsidered in detail by several teams (Nobili, Turolla & Zane, 2011; Zane et 
ah, 2011b). The main point is that single photon pair production requires photons with 
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energy higher than the threshold value, ~ 1 MeV. Therefore, in a region dominated by 
resonant scattering, pair creation occurs in two steps: (i) a seed photon with energy e ~ 1 
keV is up-scattered by a relativistic particle with 7 = 'jres ~ /^){.B/Bq) ~ 1000, 

where is the charge Lorentz factor at resonance, gaining a considerable energy 
e' ~ 7 ^g^e/(l -|- 7 rese/meC^); (ii) quite immediately, the high-energy photon converts 
to a pair, via single photon pair production. As discussed by Zane et ah (2011b), 
the pair-dominated region is very thin and located just above the star’s surface where 
B > O.OSi^Q. Here a quasi-equilibrium conhguration is reached with a pair multiplicity 
~ L/\acc,res of a fcw, where L is the length of the held line and \acc,res is the distance 
travelled by a charge before reaching a Lorentz factor '^res- Screening of the electric 
held limits the potential drop to e^o/^^eC^ ~ Ires ~ 500(i?/i?Q) and the maximum 
Lorentz factor is '^res- Charges undergo only a few scatterings with thermal photons, but 
they lose most of their kinetic energy in each collision. In practice, the result is that a 
steady situation is maintained against Compton losses because electrons and positrons 
are re-accelerated by the electric held before they can scatter again. The newborn 
charges are accelerated by the huge electric held that permeates the magnetosphere up 
to a limit value, so that a cascade of pairs is generated. This runaway process limits the 
value of 7 to the threshold value for pair production, ~1000. Since pairs with 7 ~ 
are injected from this inner region into the external region, the circuit represented by 
the held lines behaves quite diherently from a double layer, allowing the current to be 
conducted with only a small potential drop (see also Beloborodov 2011). 

A detailed investigation of charge distribution in a twisted magnetosphere, based 
on analytical considerations and corroborated by numerical tests has been recently 
presented by Beloborodov (2013a,b). This work conhrms the presence of an inner 
region with intense pair production. This region consists of two parts. The innermost 
one, where B 3> Bq, is self-organized maintaining the near critical condition of pair 
production with multiplicity M ~ 1 , and here the circuit operates as a global discharge 
(i.e. the accelerating voltage, which is screened by pairs, is distributed smoothly along 
the held line). Field lines that extend to larger distance from the star’s surface enter an 
outer corona, which extends until B ~ Bq, where both scattering and pair production 
are much more efficient and M ~ 100. Here, due to efficient radiative coupling, plasma 
and radiation organize themselves into a “locked” outhow with decreasing Lorentz factor. 
This leads to the formation of an extended equatorial zone in the outer corona, where 
the how is slowed down by the combined ehect of a large radiation drag and the onset of 
a two-stream instability with consequent strong Langmuir turbulence. The pair density 
is near annihilation balance and the charges, decelerated down to mildly relativistic 
velocities, creates an opaque layer which efficiently up-scatters the soft X-ray photons 
by distorting the surface thermal spectrum. Outside this equatorial region, and further 
away in the extended external magnetosphere, charges how at ultra-relativistic velocity 
and scattering is relatively inefficient. The charge distribution is illustrated in Fig. 6 . 
Despite this being the most complete study of magnetospheric currents presented so far, 
it contains some drastic simplifying assumptions. The pair multiplicity, for example. 
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Figure 6 Sketch of an activated magnetic loop as proposed by Beloborodov (2013b). 
Relativistic particles are injected near the star where B > Bqed- Large multiplicity 
(M ~ 100) develops in the adiabatic zone B > 10^^ G (shaded in blue). The outer 
part of the loop is in the radiative zone; here the scattered photons escape and form 
the hard X-ray spectrum. The outflow decelerates (and annihilates) at the top of the 
loop, shaded in pink; here it becomes very opaque to the thermal keV photons flowing 
from the star (from Beloborodov, 2013b, ©AAS. Reproduced with permission. A link 
to the original article via DOI is available in the electronic version). 


is a constant parameter assumed a priori. This clearly affects many of the model 
results: local screening of the electric field, velocity distribution of the two species, 
development of the two-stream instability, efficiency of the radiative drag, and ultimately 
the formation of a zone hlled with slowly moving charges. The challenging problem to 
hnd a fully consistent solution of current dynamics including the interaction with the 
radiation held remains so far unsolved. 

3.3. Soft and Hard X-ray spectral modelling 

3.3.1. Soft X-ray spectral modelling The soft X-ray band (0.3 — 10 keV) is the 
energy range in which magnetar spectra are best studied, thanks to a large amount 
of observations that have been collected in more than two decades with X-ray satellites 
such as XMM-Newton, Chandra, Swift. The observed spectra are generally htted by 
a double component model consisting of a thermal component (a blackbody, at about 
~ 0.5 keV) and a steep power law (photon index ~ 3-4) (Mereghetti, 2008). In a few 
cases (a notable example being the AXP XTE J1810-197) good spectral hts are also 
obtained with two blackbodies (Halpern & Gotthelf, 2005). The thermal component, 
which often dominates in the lowest energy band, most often has an inferred emission 
region (for the best estimated distances) much smaller than the whole surface of the 
NS. These spectral hts are of course phenomenological descriptions, but they indicate 
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that, although the emission is mostly thermal, it is more complex than a blackbody at 
a single temperature. It is also interesting to note that when one compares the average 
temperature (from the thermal luminosity) of magnetars with those of other classes of 
isolated neutron stars, there is a clear correlation between temperature and magnetic 
held (Aguilera et al., 2008) and the magnetars are systematically more luminous than 
rotation-powered neutron stars of comparable characteristic age (see a discussion in 
Mereghetti et ah, 2014). The morphology of the soft X-ray pulse prohles of magnetars 
is varied. A few sources exhibit an (almost symmetric) double-peaked light curve (e.g. 
IE 2259-1-586 and 4U 0142-1-0162; Patel et ah, 2001; Woods et ah, 2004; Rea et al., 
2007c) with pulsed fraction in the range 10 — 20%. For most other magnetars, instead, 
the pulsed component is single peaked and often the pulsed fraction is high (see e.g. 
IE 1048.1-5937, XTE J1810-197, IE 1547.0-5408, SGR 0418+5729, and SGR J1822.3- 
1606 Tam et ah, 2008; Bernardini et ah, 2009, 2011a; Halpern & Gotthelf, 2011; Dib 
et al., 2012; Rea et ah, 2013a, 2012a). As originally suggested by Marsden & White 
(2001), as a general rule sources with larger spin-down rate have smaller photon index 
in the soft X-ray band. This fact has been conhrmed with more recent data and it 
appears to be valid for both persistent and transient sources in outburst, but only for 
rotational frequencies derivatives z> > 10“^"^ s“^, and with some exceptions (including 
the transients in outbursts and the recently discovered low-B magnetars, see Sec. 4.2 
and Mereghetti et ah, 2014). The long term evolution of the power law component and 
timing properties of SGR 1806-20 indicates that the same correlation between spectral 
hardness and average spin-down rate also holds for a single source (Mereghetti et al., 
2005b). On the wake of this, other correlations between the spectral hardness and the 
timing properties have been investigated. In particular, that with the dipole strength 
Bdip oc (PP)^A appears the most robust (Kaspi & Boydstun, 2010). 

It has been widely suggested that the blackbody plus power law spectral shape 
that is observed below ~ 10 keV in magnetars’ spectra may be accounted for if 
the soft, thermal spectrum emitted by the star’s surface is distorted by resonant 
cyclotron scattering (RGS) onto the magnetospheric charges. Since electrons permeate a 
spatially extended region of the magnetosphere, where the magnetic held varies, resonant 
scattering is not expected to give rise to narrow spectral lines (corresponding to the 
successive cyclotron harmonics), but instead to lead to the formation of a hard tail 
superimposed on the seed thermal bump. This model is also in general agreement with 
the hardness-P or hardness-magnetic held correlation mentioned above: stronger and 
more twisted helds yield a larger spin down rate as well a higher magnetospheric charge 
density that in turn produces a harder spectral tail. In recent years, several teams have 
tested the resonant cyclotron scattering model quantitatively against real data in the 
soft X-ray range, using diherent approaches and under diherent approximations. The 
hrst, seminal attempts in this direction were presented by Lyutikov & Gavriil (2006) who 
developed a very simplihed one dimensional model. They assumed that seed photons 
are emitted by the NS surface with a blackbody spectrum, and propagate backward and 
forward in the radial direction. A thin, plane parallel magnetospheric slab, permeated 
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by a static, non-relativistic, warm medium at constant electron density is assumed to 
exist above the star’s surface. Magnetic Thomson scattering occurs between photons 
and the charges in the slab, and the process is computed by neglecting all effects of 
electron recoil, as are those related to the current’s bulk motion. Despite being clearly 
over-simplihed, this model has the main advantage of being semi-analytical and, when 
systematically applied to X-ray data, has proved successful in capturing at least the 
gross characteristics of the observed soft X-ray spectrum (Rea et ah, 2007a,b, 2008). 

The same model has been extended by Giiver et al. (2007), who relaxed the 
blackbody approximation for the seed surface radiation and made an attempt to include 
atmospheric effects, treating the star’s surface emission like that of a passive cooler, i.e. 
using an atmospheric code akin to those originally developed for sources with purely 
thermal emission (e.g. Zavlin et ah, 1996; Zane et ah, 2001; Ozel, 2003; Potekhin, 2014, 
and references therein). This is a quite drastic and somewhat unphysical simplihcation 
for sources like magnetars, which are characterized by strong magnetospheric activity 
leading to particle back-bombardment, heat deposition and other similar effects. Despite 
this, the model has been applied to real data in an attempt to estimate the surface 
magnetic held through data htting of the soft X-ray continuum (Giiver et ah, 2008, 
2011 ; Ozel, 2013). At present the problem appears to be still open: while it is commonly 
recognized that thermal radiation from the star’s surface is likely to be different from 
a simple blackbody, either because of local reprocessing by some sort of (non passive) 
atmosphere or because the surface itself may be in a condensed state, a self consistent 
inclusion of these effects in numerical models has not yet been carried out. 

In order to perform a more physical, 3-D treatment of the RGS problem, the most 
suitable approach is to make use of a Monte Garlo technique, which is quite easy to 
code, and, when dealing with purely scattering media at moderate optical depths, 
relatively fast. The Monte Garlo scheme allows one to follow individually a large sample 
of photons, treating probabilistically their interactions with charged particles. These 
simulations have been developed by only a few teams (Fernandez & Thompson, 2007; 
Nobili, Turolla & Zane, 2008a). The numerical codes that have been developed are 
completely general, inasmuch that in principle they can handle different 3-D geometries 
(so highly anisotropic thermal maps and magnetic helds) and different radiative models 
of surface emission. On the other hand, since our understanding of these ingredients is 
still limited, in order to minimize the number of degrees of freedom, simulations were 
computed by assuming, for simplicity, that i) the whole surface emits isotropically as 
a blackbody at a single temperature, ii) the magnetic held is a force-free, self-similar, 
twisted dipole and hi) the electron velocity distribution is assumed a priori (see Sec. 3.2). 
Besides, resonant scattering was treated in the magnetic Thomson limit, which allows 
one to account for polarization (under the two stream approximation) but it neglects 
electron recoil, limiting the applicability of the results to energies up to a few tens of keVs 
{hi' < me? j'y keV, BjBq < 10). By comparing the results from the different teams, one 
may conclude that, while the general effects induced by magnetospheric RGS on primary 
thermal photons (i.e. the formation of a “thermal-plus-power-law” spectrum) are not 



Magnetars: the physics behind observations 


24 




Figure 7 Left: Computed spectra from Monte Carlo simulations for B = G, 
kT = 0.5 keV, kTf, = 30 keV, A0 = 1 and different values of Idbuik'- 0.3 (dotted), 0.5 
(short dashed), 0.7 (dash-dotted) and 0.9 (dash-triple dotted). The solid line represents 
the seed blackbody and spectra are computed at a magnetic colatitude: = 64°. 

Right: Spectrum from a single emitting patch on the star surface. The line of sight is 
at Qs = 90° and = 20° (dotted line), 140° (dashed line) and 220° (dash-dotted line). 
These three values correspond to having the emitting patch in full view (seen nearly 
face on), partially in view and screened by the star. The solid line represents the seed 
blackbody (readapted from Nobili, Turolla & Zane, 2008a, with GUP permission). 


IE 1048-5937 SGR 1806-20 



very sensitive to the assumed particle velocity distribution, the details of the spectral 
shape do, and, as we will discuss later on, this is particularly critical for the model 
predictions in the hard X-ray band. Nevertheless, in the soft X-ray band, for several 
combinations of the parameters, the general shape of the continuum is that of a thermal 
bump and a high-energy tail (see Fig. 7), which is in agreement with what is observed 
in the XMM-Netwon and Chandra spectra (below ~ 10 keV). The spectral index of 
the high energy tail changes with the parameters and, in particular, harder spectra are 
found for increasing twist angle. This was invoked as a possible mechanism to explain 
the correlated flux-hardening long term variations in some sources (e.g. Thompson, 
Lyutikov & Kulkarni, 2002; Mereghetti et ah, 2005b; Rea et ah, 2005; Campana et ah, 
2007; Nobili, Turolla & Zane, 2008a). 

The numerical spectra computed by Nobili, Turolla & Zane (2008a) have been 
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Energy (keV) Energy (IceV) 

Figure 8 Fit of the XMM-Newton spectra of the AXP IE 1048-5937 and SGR1806-20 
with the NTZ model. The panels show a joint fit of spectra taken at three different 
epochs and the fitting has been restricted to the 1 — 10 keV range (adapted from Zane 
et ah, 2009, with permission). 


implemented in XSPEC and successfully fit the XMM-Newton spectra of most of the 
magnetar sources in quiescence (NTZ model, Zane et ah, 2009, see also Fig. 8). This 
allows one to derive of the gross characteristics of the magnetosphere in such sources 
and to obtain a better estimate of the thermal component. 

Interestingly, in the case of two sources, IE 2259-1-586 and 4U 0142-1-614 , it was not 
possible to find a satisfactory fit with the NTZ XSPEC model although these spectra 
were htted by the simplified 1-D model (Rea et ah, 2007a). As suggested in Zane 
et ah (2009), a plausible cause is that the BB peak appears to be less prominent to 
the observer because the region that emits the soft seed photons is not completely in 
view. By modelling the RCS spectra under the assumption that photons are emitted 
by a single surface patch it is immediately clear that the effects of the different viewing 
angle on the spectrum are dramatic. When the emitting patch is in full view both 
the primary, soft photons and those which undergo repeated resonant scattering reach 
the observer, and the spectrum is qualitatively similar to those presented earlier on, 
with a thermal component and an extended power-law-like tail. If on the other hand 
the emitting region is not directly visible, no contribution from the primary blackbody 
photons is present (see Fig. 7, right panel). The spectrum, which is made up only 
by those photons which after scattering propagate “backwards”, is depressed and has 
a much more distinct non-thermal shape, much more similar to the one observed in 
IE 2259+586 and 4U 0142+614 . 

Unfortunately, in a realistic situation the thermal surface map is expected to be 
complex and it cannot always be reconstructed by htting the X-ray spectrum alone. As 
discussed in Albano et al. (2010, see also Bernardini et ah 2011a), a better strategy 
consists of performing a simultaneous fitting of the energy-resolved X-ray lightcurves, 
since they carry a much more defined imprint of the surface thermal distribution (see 
Sec. 4 and Fig. 18 therein). While this is not possible for all sources, transient AXPs, for 
which a set of observations spread over few years and at different flux levels are available, 
provide a spectacular laboratory for this exercise. Albano et al. (2010) were the hrst (and 
so far the only) team to present a comprehensive study of the outburst decay of the two 
transient AXPs (TAXPs) XTE J1810-197 and CXOU J164710.2-455216, reproducing 
both the spectral and pulse profile evolution based on hts with three-dimensional Monte 
Carlo simulations. This allowed them to prove the presence of distinct temperatures 
zones (up to three) at the star’s surface, some of them possibly heated by the energy 
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released during the outburst, to model their evolution during the outburst decay, and 
to constrain the viewing geometry of the sources, i.e., the inclination of the line of sight 
and the magnetic axis with respect to the rotation axis. 

A similar conclusion concerning the need to investigate the pulse prohle behaviour 
when trying to reconstruct inhomogeneous surface thermal maps was reached by Perna 
et al. (2013). These authors considered the case in which the temperature anisotropy is 
observed also in quiescence, as is expected if complex magnetic held components in the 
NS crust and interior make heat transport from the core outward highly anisotropic. 
They used state-of-the art numerical codes (Viganb et al., 2012, 2013) for the coupled 
magneto-thermal evolution of neutron stars and computed the expected pulse prohles 
and spectra (under the assumption of blackbody emission) for a range of magnetic 
conhgurations. Particularly compelling is the hnding that, while in presence of purely 
dipolar helds the pulse prohle is always double-peaked and with a relatively low pulsed 
fraction, when strong toroidal components are present the pulse fraction can exceed 50- 
60% and the pulse prohle can be single peaked (as often observed in AXPs and SGRs). 
Moreover, if the simulated spectra are htted with a highly absorbed BB model, only 
relatively concentrated hot peaks are visible, so that the inferred BB radius turns ont 
to be much smaller than the NS radius (even as low as 1-2 km, see Fig. 9). Strong 
toroidal crustal B-held components, coupled with large absorption column densities 
(> 10^^ cm“^, see e.g. Esposito et al. 2008), can therefore explain the small caps 
very often reqnired by the spectral hts of AXPs and SGRs. Even smaller (snb-km) 
hot spots are measured in certain sources (e.g. GXO J164710.2-455216, see e.g. Israel 
et al. 2007), but they look more likely to be produced by particle bombardment and 
heat deposition from currents highly concentrated in twisted magnetic polar bnndles 
(Beloborodov, 2009; Turolla et al., 2011) which emerge from the crust, rather than 
anisotropic internal heat transfer. In this respect it is worth mentioning the work 
by (Bucciantini et ah, 2015), who presented a comprehensive and detailed parameter 
study, in general relativity, of the role that the cnrrent distribution investigating several 
eqnilibrium global held conhgurations derived using a Grad-Shafranov approach (Pili et 
al., 2015). These anthours fonnd that the structnre and strength of the magnetic held at 
the snrface is strongly inhnenced by the location and distribntion of currents inside the 
star, with the result that the surface held can easily be dominated by higher multipoles 
than the dipole. This means that in some cases the magnetic held at the equator can 
be even mnch higher or mnch smaller than the value of the held at the pole and implies 
that signatures observed in features originating at or near the surface might diher from 
the expectations of a dipole dominated model, while observations of processes related 
to the large scale held, as spin-down, will not (see also the discussion in Sec. 4.2). 

3.3.2. Hard X-ray spectral modelling Hard X-ray observations with the INTEGRAL., 
RXTE and Suzaku satellites have shown that in some magnetar candidates (namely 4U 
0142+614, IRXS J1708-4009, IE 1841-045, IE 2259+586, SGR 1806-20, SGR 1900+14; 
Kniper et al., 2004, 2006; Mereghetti et al., 2005a; Molkov et ah, 2005; Gotz et al., 2006; 
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Figure 9 Temperature surface distribution for sources with a magnetic field at birth 
which has both a poloidal and a toroidal component. For ages typical of magnetars 
(10^ — 10^ yrs) the expected thermal map consists of a hot polar spot and an extended 
colder region; the latter may be undetected if the source is highly absorbed (from Perna 
et al., 2013, with OUP permission). 


Enoto et ah, 2011) a large fraction of the total quiescent flux is emitted at energies well 
above ~ 20 keV. These “hard” tails have a non-thermal (PL) character, and extend 
up to a few hundreds of keV. Pulsed phase spectroscopy has been performed for a few 
sources, although with limited statistics due to the low number of counts, revealing that 
the emission is likely characterized by different components that emerge at different, 
and sometimes only in limited, phase intervals (den Hartog et al., 2008a,b). 

This discovery came quite unexpectedly and suggests that a new magnetar 
characteristics may be that a considerable fraction of their bolometric luminosity is 
emitted in the hard, rather then in the soft X-rays. In fact, limits on the non-detected 
sources are not deep enough to exclude the presence of a similar hard tail. The hrst 
observations of emission at > 20 keV revealed a difference between the (at that time 
separated) AXPs and SGRs (Gotz et ah, 2006): in the Integral and RXTE data the 
AXPs’ hard power law is considerably harder than that observed below 10 keV, while 
SGRs’ hard spectra are considerably steeper. More recent observations with NuStar 
indicate that the division is not as sharp and magnetar sources show a varied behaviour, 
which is consistent with the fact that SGRs and AXPs are now considered a single class 
(An et ah, 2014b). NuStar, which has a sensitivity roughly two orders of magnitude 
better than previous missions in this energy band and a high angular resolution, is 
currently observing a selected sample of magnetars as part of its priority A targets (see 
An et ah, 2014b, for a review and references therein). One of the major goal is to study 
the spectral location of the soft/hard X-ray turnover, which is expected to correlate with 
Bdip (Kaspi & Boydstun, 2010). Interestingly, while for some sources the new NuStar 
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Figure 10 Left: XMM-Newton and INTEGRAL spectra of magnetars (from Gotz et al., 
2006, reproduced with permission (c)ESO). Note the different behavior of SGRs (two 
top panels) and AXPs: in the latter sources the spectra turn upward above 10 keV, 
while in the SGRs the spectra steepen. Right: Broad band spectral energy distribution 
of 4U 0142+614 (from Kuiper et ah, 2006, (c)AAS. Reproduced with permission. A link 
to the original article via DOI is available in the electronic version). Both the total and 
the pulsed emission are indicated. Both hgures are readapted from Mereghetti (2008), 
with kind permission from Springer Science and Business Media. 


results are in agreement with the Integral and RXTE ones, in some other cases (e.g. 
IE 2259+586) it appears evident that the hard band spectrum for the total emission is 
not as hard as the pulsed one. Despite that, the analysis of IE 2259+586 data conhrmed 
that an additional component, such as a power law, is needed to describe the emission 
in the hard X-ray band (Vogel et ah, 2014). This suggests that at least in some sources 
the non thermal mechanisms responsible for the emission in the soft and hard X-rays 
are distinct, or that the charge properties in the two emission regions are different. 

Observations at higher energy with Gomptel and Fermi LAT failed to detect 
magnetar emission, implying the presence of a spectral break above a few hundred 
keV (Kuiper et ah, 2006; den Hartog et ah, 2006; Sasmaz Mus & Gogu§, 2010). The 
only exception reported so far is a possible Fermi Large Area Telescope (LAT) detection 
of 7 -ray pulsations above 200 MeV from the AXP IE 2259+586 (Wu et ah, 2013) which 
however still needs a robust conhrmation. Some example of few magnetar spectral 
energy distributions in the soft/hard X-rays are shown in Fig. 10. 

The mechanism responsible for the high energy emission is still poorly understood. 
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at least in its quantitative details. High energy emission from currents moving in the 
highly magnetized magnetosphere is expected for a number of reasons. Thompson 
& Beloborodov (2005) originally suggested that hard X-rays may be produced either 
by thermal bremsstrahlung in the surface layers heated by returning currents, or by 
synchrotron emission from pairs created higher up (~ 100 km) in the magnetosphere. 
A further possibility, according to which the soft y-rays may originate from resonant 
up-scattering of seed photons on a population of highly relativistic electrons, has been 
proposed by Baring & Harding (2005, 2008). As mentioned earlier, most of the RCS 
models computed with Monte Carlo simulations are based on the assumption that 
scattering can be treated in the Thomson approximation, which limits their validity 
up to a few tens of keV. Instead, a proper investigation of the effects of electron recoil 
and of multiple scatterings from high energy photons demand the use of the full QED 
cross section. Relatively simple expressions of the QED cross section at resonance, in a 
form that is simple to include in Monte Carlo simulations, have been computed by Nobili, 
Turolla & Zane (2008b) and a few examples of the emerging spectra, computed under 
the assumption of self-similarly distributed twist and constant electron velocity, have 
been discussed in Zane et ah (2011a). These simulations show that, if magnetospheric 
electrons are mildly relativistic, when considering self-consistently electron recoil and 
QED effects the spectrum exhibits a break at a few hundred of keV. This is due to the 
fact that, in order to populate the hard energy tail, soft seed photons need to experience 
a series of successive scatterings, each characterized by a limited energy gain because 
the Lorentz factor of electrons is only 7 ~ a few. In parallel, the efficiency of the 
QED cross section decreases with increasing energy (or, since the process is resonant, 
with increasing magnetic held) and the combination of these two effects leads to the 
appearance of the spectral break. The energy of the break depends on the effect of the 
cumulative scatterings and is sensitive to the details of the magnetic held topology and of 
the currents’ distribution and therefore cannot be predicted a priori nor estimated using 
a simple expression. On the other hand, in the case in which magnetospheric electrons 
are ultra-relativistic, the energy gain per scattering is so large that the hard tail becomes 
efficiently populated after just a few scatterings. In this case, now independently of the 
details of the cross section and magnetic topology, the spectral tail is predicted to be 
quite hat and unbroken, even up to > 1000 keV (see Fig. 11). 

Although these studies are extremely useful in shedding light on the basic behaviour 
of the QED scattering process, the assumptions of constant charge velocity and self¬ 
similar magnetic twists are clearly two major oversimplihcations. This is crucial when 
trying to mimic the hard X-ray emission, since the responsible emitting region is likely 
to constitute quite a large portion of the whole magnetosphere. We already mentioned 
that deep INTEGRAL observations of two AXPs IRXS J1708-4009 and 4U 0142-1-61 
have revealed several diherent pulse components (at least three) in the hard X-rays, 
with genuinely diherent spectra (den Hartog et ah, 2008a,b) and a quite spectacular 
phase-dependence. The hard X-ray spectrum gradually changes with phase from a soft 
to a hard power law, the latter being signihcantly detected over a phase interval covering 
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Figure 11 Monte Carlo spectra for B = 10^^ G and A0 = 1 computed by using the 
full QED expression of resonant scattering. Left: mildly relativistic electrons (7 = 1.7). 
Right: highly relativistic electrons (7 = 22) (reprinted from Zane et ah, 2011a, Copyright 
2011, with permission from Elsevier). 


~l/3, or more, of the spin period. This richness in the phenomenology requires more 
complex magnetospheric topologies to be explained. Unfortunately, despite many efforts 
having been devoted to the development of techniques for solving the force-free equation, 
V X B = a{x)B, no general, affordable method has been presented so far. Pavan et ah 
(2009) developed a general method to generate twisted, higher-order multipoles solving 
the Grad-Shafranov equation, and analyzed in detail quadrupolar and octupolar helds. 
The case of an octupolar field has a special interest because it can be used to mimic 
a twist localized in a region close to the magnetic pole(s), and hence to investigate 
the consequence in the expected spectra. Model (Monte Carlo) spectra and lightcurves 
have been presented for the cases in which the twist is conhned to one or both polar 
regions (each region has semi-aperture of ~ 60deg), by assuming that only the polar 
lobes have a non-vanishing shear while the equatorial belt is potential. Interestingly, 
a conhguration with a twist confined to a single lobe has been found to be capable of 
qualitatively reproducing the main features of the high-energy emission observed with 
INTEGRAL from the AXPs IRXS J1708-4009 and 4U 0142-1-61, in particular the large 
variation in the pulsed fraction at different energy bands, and a hard tail which is quite 
pronounced at the peak of the pulse but depressed by almost an order of magnitude at 
pulse phases close to the minimum of the hard X-ray lightcurve. 

More recently, an alternative to the (mathematically simple) self-similar models, 
has been presented by Vigano et ah (2012), who discussed the effects of more 
realistic magnetic field geometries on the synthetic Monte Carlo spectra. They 
presented a numerical method to build general force-free held magnetic conhgurations, 
starting from an arbitrary, non-force-free poloidal plus toroidal held and employing 
artihcial dissipation to remove the non-parallel currents. In particular, they considered 
conhgurations in which the currents are concentrated in a bundle along the polar axis, 
as expected for a spatially-limited twist. In this case the pulse fraction is larger with 
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Figure 12 Left and center panel: magnetic field and current distribution for a high 
helicity model with current and twist concentrated in a closed bundle near the equatorial 
region and a j-bundle near the southern semi-axis. This conhguration is likely more 
realistic than the self-similar models, since currents are more concentrated near the 
axis. The two panels show \B\ (grey logarithmic scale) with superimposed the scattering 
surfaces for photons of 1, 3 and 5 keV and the current intensity | J|/c (gray linear scale) 
in units of 10^^ G/r^s- Right: Corresponding synthetic spectra computed with a Monte 
Carlo simulation. Different curves correspond to four different viewing angles; the seed 
blackbody is shown for comparison as a solid line, (readapted from Vigand et ah, 2012, 
with permission; a link via DOI to the original version is available in the electronic 
version of this paper). 


respect to that of self-similar models, and the spectrum observed at different angles 
varies in a much more irregular way (see Fig. 12). Instead of a simple PL, it shows 
different spectral bumps the relative importance of which can vary by one order of 
magnitude or more at different colatitudes. The different spectral components inferred 
from the data may therefore be due to these bumps. Even if a real ht has not been 
attempted, we may speculate that this is qualitatively in line with the observations of 
the AXPs IRXS J1708-4009 and 4U 0142+61. 

As mentioned earlier, a further poorly known ingredient of all these simulations is 
the charge velocity distribution and probably the most detailed solutions published have 
been presented by Beloborodov (2013a,b). In this scenario, the electron-positron flow 
decelerates as it propagates away from the neutron star surface (due to Compton drag 
in the resonant scattering region), then it reaches the top of the magnetic loop where it 
annihilates. While computed with a Monte Carlo simulation, the corresponding spectra 
show a distinct peak at E > 1 MeV (so far unobserved), the position of which is however 
strongly dependent on the viewing angle (or, equivalently, on the magnetic colatitude 
at which the spectrum is emitted, see Fig. 13). 

A well dehned change in the power law slope from soft to hard X-rays is seen, which 
makes these spectral models incapable of explaining observations of sources that have 
a similar slope below and above ~ 10 keV (historically these were referred as SGR-like 
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Figure 13 Left: Spectrum emergent from scattering onto charges populating an 
axisymmetric j-bundle. Different lines represent the spectrum as observed at four 
different angles with respect to the magnetic axis: 15°, 30°, 60°, and 90°. The spectrum 
is normalized to the total luminosity of the j-bundle. (from Beloborodov, 2013b, ©AAS. 
Reproduced with permission. A link to the original article via DOI is available in the 
electronic version). 


spectra), but instead in good qualitative agreement with the observed spectra in which 
a marked turnover is present (historically, AXP-like spectra). In an attempt to perform 
a quantitative comparison without to resort to the time-consuming computation of a 
complete table of Monte Carlo models, Hascoet et ah (2014) developed a simplihed model 
by assuming the same charge velocity distribution but by computing the spectra through 
a simple calculation of the angle-dependent emissivity alone. The magnetosphere is 
assumed to be axysymmetric and dipolar, apart from the presence of a current-carrying 
region (the j-bundle), which is hlled with the electron-positron flow. These spectra have 
been then applied to multi-year hard X-ray observations of 4U 0142-1-61, IRXS J1708- 
4009 and IE 1841-045, Ending that they successfully reproduce the emission observed 
above ~ 10 keV in both the phase-average and phase-resolved spectra. Unfortunately, 
the authors found that the model predictions cannot be self consistently applied to 
explain also the data below 10 keV (probably because of the lack of the particle back- 
bombardment efiects in the simulation). These predictions have therefore been ignored 
in the fits and instead a blackbody component (s) was used to account for the soft X- 
ray part of the spectrum. A concern regarding these fits is the fact that the spectra 
used by Hascoet et ah (2014) do not include photon splitting, which is the main factor 
responsible for the bump visible in the Monte Carlo spectra at ~ 1 MeV (Beloborodov, 
2013b). Although the simplified and complete calculations are in good agreement below 
~ 400 keV, they are expected to be markedly different at higher energies. Actually, 
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there is no guarantee that, should the spectra be computed self-consistently with the 
Monte Carlo code, the high energy CGRO-Comptel upper limits are not violated (either 
always or for the derived geometry). This casts some doubts on the robustness of the 
study and on the constraints on the viewing geometry reported by Hascoet et al. (2014) 
htting only the hard X-ray part of the spectrum. 

It is interesting to note that spectra presented by Wadiasingh et al. (2013), using 
a new Sokolov and Ternov formulation of the QED Compton scattering cross section in 
strong magnetic helds and accounting for spin-dependent effects at resonance, also show 
a spectral cut-off whose energy is critically dependent on the observer viewing angles 
and electron velocity, with substantial emission expected up to 1 MeV except for very 
selected viewing angles. 

At the present, unfortunately, what causes the high energy emission and its richness 
in phenomenology is still an open issue: models have been computed, but they depend 
dramatically on a large number of degrees of freedom for the magnetospheric setting at 
large scale. New and future missions such as Astro-H, NuSTAR, and possibly LOFT 
would provide the possibility of collecting simultaneous soft and hard spectra, to study 
the correlation between the variability in the two bands, to perform high resolution 
pulsed phase spectroscopy and to reveal the detail of the slope turn-over (when present) 
in the soft-to-hard emission. In turn, this will provide a powerful tool to break the 
model degeneracy and an unprecedented insight into the details of the held and current 
distributions. 

3.4- X-ray polarization 

As discussed in Sec. 3.3.1, comparison of RCS models with the soft X-ray spectral data 
of magnetar candidates has proved quite successful. However, spectroscopy alone cannot 
provide complete information on the physical properties of the magnetosphere, due to 
the inherent degeneracy in the RCS model parameters. Moreover, computed spectra 
are rather insensitive to the source geometry, although in principle they do depend on 
the angles that the line of sight and the magnetic axis make with the star’s rotation 
axis (Nobili, Turolla & Zane, 2008a; Zane et al., 2009). Although a simultaneous ht of 
both the (phase-averaged) spectrum and the pulse prohle is effective in this respect, at 
least for TAXPs (Perna & Gotthelf, 2008; Albano et ah, 2010; Bernardini et ah, 2011a), 
polarization measurements at X-ray energies would provide an entirely new approach 
to the determination of the physical parameters in magnetar magnetospheres. 

X-ray radiation from a magnetar is expected to be polarized for essentially three 
reasons: i) primary, thermal photons, coming from the star’s surface, can be intrinsically 
polarized, because emission favors one of the modes with respect to the other; ii) 
scattering can switch the photon polarization state; and iii) once the scattering depth 
drops, the polarization vector changes as the photons travel in the magnetosphere (the 
so called “vacuum polarization”, Heyl & Shaviv, 2000, 2002, see also Harding & Lai 
2006). 
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Figure 14 Contour plots for number of counts (in arbitrary units; left column), 
polarization fraction (middle column) and polarization angle (right column) as functions 
of the photon energy and cos 6* for different values of the twist angle and the electron 
bulk velocity: A0 = 1.3 rad, (3 = 0.3 (top row) and A0 = 0.7 rad, (3 = 0.4 (bottom row). 
In both cases it is i?p = 5 x 10^^ G (from Taverna et ah, 2014, with OUP permission). 

Although a preliminary analysis was already contained in Fernandez & Thompson 
(2007) and Nobili, Turolla & Zane (2008a), a detailed study of the polarization properties 
of magnetar radiation in the X-ray band has been presented by Fernandez & Davis (2011) 
and, more recently, by Taverna et ah (2014), by means of Monte Carlo simulations. 
Phase-averaged as well as phase-resolved results indicate that the linear polarization 
fraction IIl and the polarization angle Xpoi, are very sensitive to the magnetospheric 
twist angle A0 and the charge velocity (3, and also to the geometric angles x and 
This allows one to remove the A(j)-(3 degeneracy which spectral measures alone cannot 
disambiguate. An example is shown in Fig. 14 where the photon spectrum, polarization 
fraction and polarization angle are plotted as functions of viewing angle and energy 
for two choices of the model parameters. While the photon spectrum is practically 
the same, the pattern of 11/, and especially Xpoi is markedly different in the two cases. 
According to the simulations by Taverna et ah (2014), polarimetric measurements in 
bright magnetar candidates, like the AXP IRXS J1708-4009, are within reach of recently 
proposed polarimeters which should hopefully fly in future missions. 

3.5. IR/optical emission 

As mentioned in Sec. 1, variable IR counterparts have now been identified, or in some 
cases proposed, for a number of magnetars (Mereghetti, 2011; Israel & Rea, 2014). In 
addition, three magnetars have been detected in the optical band, the two AXPs 4U 
0142-1-61 (HuIIeman et ah, 2000) and IE 1048.1-5937 (Durant & van Kerkwijk, 2005), 
and the SGR 0501-1-4516 (Fatkhullin et ah, 2008; Dhillon et ah, 2011). The origin of this 
emission is still under debate, and the main dispute is whether it is due to the presence 
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of a fossil disk (Perna et al., 2000) or if it has a magnetospheric origin (Eichler, Gedalin 
& Lyubarsky, 2002; Beloborodov & Thompson, 2007). 

AXP 4U 0142+61 is the only persistent magnetar for which the optical/near IR 
(NIR) spectrum was measured, although the faintness of the source required extensive 
observational data (Hulleman et ah, 2000). In the case of this source, an IR “excess“ 
(or “flattening “) was clearly detected with respect to the extrapolation of the additional 
black body used to account for the optical emission (Israel et ah, 2004, 2005a). This may 
suggest that the IR emission is due to a distinct spectral component, a conclusion that 
might well hold also for other AXPs, considering the similarity of their IR magnitudes 
and Fx/Fi^ ratios (Durant et ah, 2011). 

The optical/NIR spectrum of 4U 0142+614 was found to be well htted by a multi¬ 
temperature (700-1,200 K) thermal model, leading to the suggestion it originates in an 
extended disk or shell (Wang et ah, 2006). If this detection of a fallback disk is real, 
it would be the the first direct evidence for supernova fallback in any context. On the 
other hand, disks may then be ubiquitous while, despite intense campaigns, no similar 
direct evidence has been found for other sources (see also Posselt et ah, 2014, for a 
report on recent deep limits on fallback disks). The only indirect evidence is one source, 
in which there is a detected correlation between the NIR and X-ray fluxes (Tam et ah, 
2004), which, as pointed by Wang et al. (2006), may indicate that the IR emission arises 
in an X-ray-heated debris disk. 

On the other hand, very deep optical and NIR observations of the held of the low-B 
magnetar SGR 0418+5729, which is the nearest and least extincted magnetar known, 
failed to detect the source counterpart (Durant et ah, 2011). Shallower observation of 
the held of SGR 0418+5729 with the new Gran Telescopio Ganarias 10.4-m telescope 
and with the William Herschel Telescope where also taken closer to the onset of the 
outburst (Esposito et ah, 2010; Rea et ah, 2013a) and only gave upper limits on the 
countarpart. This negative result is in better agreement with a magnetospheric origin 
interpretation, in which case the IR/optical hux is expected to be fainter for lower held 
strengths. Moreover, a magnetospheric scenario is more likely to explain the fact that, 
when detected, the observed optical emission is pulsed at the pulsar period, with pulsed 
fractions > 50%, i.e., higher than in soft X-rays (Kern & Martin, 2002; Dhillon et ah, 
2009, 2011). 

Magnetospheric emission in the IR/optical is expected from the inner region of 
the magnetosphere (Beloborodov & Thompson, 2007; Zane et ah, 2011b). As discussed 
in Sec. 3.2, this region is expected to be pair-dominated, and the Lorentz factor of 
the electrons and positrons is likely to be frozen at the threshold for pair production 
(7 ~ 500 — 1000). In order to estimate the amount of curvature radiation, Zane et 
al. ( 2011 b) developed a simple geometrical model accounting for misalignment between 
the observer line of sight and the magnetic axis. Although the model is too simple to 
allow for a proper spectral fitting, it is detailed enough to make a prediction about the 
energetics and the computation indicates that curvature radiation is sufficient to explain 
the amount of observed IR/optical flux, at least if a particle bunching mechanism is 
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Figure 15 Model spectra for different values of the viewing angle and B = 4.3 x 10^^ G. 
The XMM-Newton X-ray spectrum of IRXS J1708-4009 is from Rea et ah (2008, red solid 
line). The AXPs IR/optical data are from Durant & van Kerkwijk (2005, 4U0142-1-614 
and IE 1048-5937) and Mignani et ah (2007, XTEJ1810-197 and IE 2259-1-586). The 
adopted distances for de-reddening are 5 kpc (4U0142-1-614, IRXS J1708-4009), 3 kpc 
(IE 1048-5937, IE2259+586), 4 kpc (XTE J1810-197), 8.5 kpc (IE 1841-045). Curvature 
emission spectra have been computed accounting for particle bunching (Fig. 2 from Zane 
et ah, 2011b, With kind permission from Springer Science and Business Media) . 


efficient (see also Beloborodov & Thompson 2007 and see Fig 15). This is not unlikely, 
since many models have been suggested to explain the origin of the interactions which 
push particles together and can lead to the formation of bunches of charged particles 
localized in phase-space. The most promising explanation seems to be connected with 
plasma instabilities, like the two-stream (electron-positron/electron-ion) instability. We 
notice that this mechanism does not affect the emission in other bands. In fact, in order 
the process to be effective, it has to have N = nj^ > 1, where Ib is the size of the bunch 
and rie is the electron density, which in turns means that only emission at frequencies 
< i^co = c is efficiently amplihed. Also, a low-energy cut-off is present because of 
the strong absorption below the electron plasma frequency. 

Israel et ah (2005b) proposed a link between the IR and hard X-ray spectrum 
of AXPs, and correspondingly between AXPs/SGRs and radio-pulsars, based on the 
analysis and comparison of their broad band energy spectra. These authors pointed out 
that, similar to the case of a number of young radio pulsars such as the Vela, it is possible 
to bridge the IR to y-ray emission of AXPs with a power-law with index of about 0.5-0.6 
(see Fig. 16). This peculiar similarity, for classes of neutron stars with otherwise quite 
different emission properties, may naturally explain the IR excess or flattening in the 
spectra of AXPs and points toward a similar origin for the less energetic bands and the 
hard X-ray emission. This scenario may be unambiguously proven through simultaneous 
studies of the correlated variability in the two bands, which have unfortunately so far 
been hampered by the lack of gamma-ray observatories sensitive enough to these faint 


sources. 
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Figure 16 Broad band energy spectrum of SGR 1806-20 (triangles), the AXP IE 
1841-1-045 (circles) and the radio pulsar Vela (squares). In the case of SGR 1806-20 and 
IE 1841-1-045 high energy data are taken from Mereghetti et ah (2005a,b). Absorbed 
and unabsorbed IR fluxes (Ay = 29±2, 5th October 2004 NAGO observation) are shown 
in the case of SGR 1806-20, unabsorbed (Ay = 13 ± 1) IR fluxes are instead reported 
for the likely candidate of IE 1841-1-045 (circles; Testa et ah, 2008). All the data for 
Vela are taken from Kaspi (2006). Solid curves (continuous, stepped and dot-stepped) 
are the unabsorbed fluxes, for the blackbody plus power-law model used to £t the high 
energy part of the spectra (courtesy G.L. Israel, from Israel et ah, 2005b, reproduced 
with permission (c)ESO). 

4. Transient Magnetars 

The persistent X-ray emission of a number of SGRs/AXPs has been known to be variable 
all along, with typical flux variations of a factor of a few over a timescale of days to 
months, often in coincidence with periods of enhanced bursting activity (see Sec. 1; Rea 
& Esposito, 2011). The hrst evidence that the luminosity of SGRs/AXPs can change 
much more dramatically came from observations of SGR 1627-41 (Woods et ah, 1999b) 
and AX J1845.0-0300 (Torii et ah, 1998; Vasisht et ah, 2000, although the latter source 
is only a candidate magnetar with no detection of P). 

It was not until 2002, however, that the existence of a new class of magnetar sources 
with much more extreme variability was realized, thanks to the discovery of the hrst 
transient AXP, XTE J1810-197 (Ibrahim et ah, 2004). At present eleven transients are 
known and, remarkably, they almost make up all of the new magnetars observed in the 
last 10 yrs (the exception being GXOU J171405.7-381031, Halpern & Gotthelf, 2010). 
The main properties of transient magnetars are listed in table 1^. 

+ Only sources for which the peak luminosity is > 10 times the quiescent one have been included. 
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Table 1 Transient magnetar sources^ 


Source 

P 

(s) 

o 

1 

B 

(1014 G) 

D 

(kpc) 

CXOU J1647-4552 

10.61 

< 0.04 

< 0.7 

3.9 

XTE J1810-197 

5.54 

0.77 

2.1 

3.5 

SGR 0501-^4516 

5.76 

0.59 

1.9 

52.0 

SGR 0418-^5729 

9.08 

0.0004 

0.06 

2.0 

SGR 1833-0832 

7.56 

0.35 

1.6 

- 

PSR 1622-4950 

4.33 

1.7 

2.2 

9.0 

IE 1547-5408 

2.07 

4.77 

3.2 

4.5 

Swift J1822.3-1606 

8.4 

0.02 

0.14 

1.6 

SGR 1627-41 

2.59 

1.9 

2.2 

11.0 

Swift J1834.9-0846 

2.48 

0.80 

1.4 

4.2 

SGR J1745-2900 

3.76 

1.38 

2.3 

8.3 


These sources are characterized by a sudden hrs) increase of the X-ray flux, by a 
factor 10-1000 over the quiescent level, accompanied by the emission of short bursts. 
This active phase, commonly referred to as an outburst, typically lasts ~ 1 yr, during 
which the flux declines, the spectrum softens and the pulse prohle simplihes. Fig. 17 
shows the decay of the X-ray flux for several transient magnetars. Some objects have 
undergone repeated outbursts (SGR 1627-41, IE 1547-5408) and the decay pattern 
is often different from source to source (and even between outbursts from the same 
source). Outbursts, besides revealing new magnetars which in quiescence are too faint 
to be detectable, or which passed unnoticed among the host of unclassihed, weak X- 
ray sources, have also occurred in persistent sources like IE 2259-1-586 or IE 1048.1- 
5937, although drawing a precise line between outbursts and less extreme variability is 
somewhat haphazard. The group of transient sources also harbours the two peculiar 
“low-£eld“ magnetars, SGR 0418-1-572 and Swift J1822.3-1606 (see Sec. 4.2), and the 
recently discovered source in the Galactic Gentre (see Sec. 4.3). 

4-1- Outburst Models 

A common feature of all observed outbursts is the presence in the X-ray spectrum of 
one (or two) thermal component(s) at higher temperature (~ 0.3-0.9 keV) with respect 
to that associated with the cooling star surface during quiescence (~ 0.1-0.2 keV). The 
(radiation) radius of these hotter regions is fairly small (< 1 km) and usually decreases 
in time as the outburst subsides, when the temperature also declines (e.g. Albano et 
al. 2010, Rodriguez et al. 2014 for XTE J1810-197, GXOU J1647-4552; Rea et ah 2009, 
Gamero et al. 2014 for SGR 0501-1-4516; Rea et al. 2013a for SGR 0418-1-572; Rea et 
al. 2012a for Swift J1822.3-1606; Israel et al. 2010 for 1E1547.0-5408, and references 
therein). The variation in time of the spectrum, pulse profile and size of the emitting 

+ Data from the McGill magnetar catalogue (Olausen & Kaspi, 2014). 
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Time (days from burst activation) 


Figure 17 Flux evolution over the first ~200 days of all magnetar outbursts (only if 
observed with imaging instruments, and for which this period span is well monitored). 
Fluxes are reported in the 1-10 keV energy range, and the reported times are calculated 
in days from the detection of the hrst burst in each source. See Rea & Esposito (2011) for 
the reference of each reported outburst. (From Rea, 2014c, Copyright ? 2014 WILEY- 
VCH Verlag GmbH & Co. KGaA, Weinheim. Reproduced with permission.) 


regions for the AXP XTE J1810-197 during the the outburst decay is shown in Fig. 18. 

This has been interpreted as due to some form of heat deposition in a limited 
region of the star surface which then cools and shrinks. Until now, however, the heating 
mechanism has not been unambiguously identihed. One possibility is that energy is 
injected deep in the crust, e.g. because of magnetic dissipation, and then flows to the 
surface, as hrst suggested by Lyubarsky, Eichler & Thompson (2002). 

Pons & Rea (2012) developed a quantitative model for the outburst evolution by 
simulating the thermal relaxation of the neutron star in response to an impulsive energy 
injection in the star crust. They found that most of the energy is released in the form of 
neutrinos, unless injection occurs in the outer crust (p < 3 x 10^^ g cm“^). The successive 
evolution depends mostly on the energy input in the outer crust, Eqc- However it 
has to be Eqc ^ 10^°ergs“^ to produce any visible effect on the surface and, at the 
same time, Eqc is bounded from above at ~ 10^^ergs“^ because any excess energy is 
efficiently radiated away by neutrinos produced in the heated crust. This limits the 
surface temperature to ~ 0.5 keV. Because heat transport occurs mostly along the held, 
which is predominantly radial in the outer layers, the size of the hot spot on the star 
surface remains nearly constant in time (which may be problematic in explaining the 
observed shrinking of the heated region). 
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Figure 18 The evolution of the spectrum, pulse profile (top panels) and size of the 
emitting regions (expressed as fraction of the entire surface, bottom panel) during the 
outburst of XTE J1810-179. Results are for the three temperature model (hot cap, warm 
corona and cool surface) of Albano et ah (2010); the cold temperature is fixed at 0.15 
keV. The adopted spectral model is the superposition of two NTZ’s (see Sec. 3.3.1) at 
the hot and warm temperature (adapted from Albano et ah, 2010, ©AAS. Reproduced 
with permission. A link to the original article via DOI is available in the electronic 
version). 


Results were successfully applied by Rea et al. (2012a) to fit the outburst decay in 
Swift J1822.3-1606 over the entire period covered by their observations, ~ 250 d after 
the hrst burst that led to the discovery of the source. The case of SGR 0418-1-5729, 
for which a much longer time coverage is available (~ 1200 d), is, however, much less 
conclusive in this respect (Rea et ah, 2013a). The calculated flux in the 0.5-10 keV 
band systematically overestimates the observed one at later times (> 400 d), when 
the luminosity suddenly drops and the hotter blackbody (initially at kT ~ 0.9 keV) 
disappears leaving only a cooler component at ~ 0.3 keV. A similar effect has been 
recently reported in the decay curve of CXOU J1647-4552 around 1000 d (Rodriguez et 
ah, 2014) and in Swift J1834.9-0846 (Esposito et ah, 2013). 

Alternatively, heating of the surface layers may be produced by currents flowing in 
a twisted magnetosphere as they hit the star. As discussed in §3.1, once implanted by 
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crustal displacements, a twist must necessarily decay in order to supply the potential 
drop required to accelerate the conduction current (Beloborodov & Thompson, 2007). 
The evolution of an untwisting magnetosphere proceeds through the expansion of a 
potential region, where V x B = 0, which progressively conhnes the twist to a limited 
bundle of current-carrying held lines (the j-bundle), until the twist is completely erased 
(Beloborodov, 2009). Fractures, or plastic deformations, most likely affect only a limited 
area of the star crust, so the twist is expected to involve only the bundle of held lines 
whose footpoints are anchored in the displaced region. As the magnetosphere untwists, 
the area covered by the j-bundle shrinks and the luminosity decreases. The rate of 
Ohmic dissipation, which depends on the initial twist conhguration and on how the 
twist angle A(j) evolves in time {A(j){t), is not necessary monotonic, the twist may hrst 
decrease and then increase to a maximum value ~ 1 rad). The simplest model gives 

^ “ ““ (lo^) (lo^) (to^) (10) 

where V is the discharge voltage and 6j-b is the angular extension of the j-bundle 
on the star’s surface (this is valid for a small polar bundle with maximal twist, see 
Beloborodov, 2009). Since heat is unlikely to leak outside the cap at the base of the 
j-bundle, the area of the X-ray emitting region is ~ sin^ ^j-b- A quite dehnite 
prediction of the model is, then, that the luminosity is proportional to the square of the 
emitting area throughout outburst evolution (this still holds, at least approximately, 
also for more elaborate versions). A spatially-limited twist can also explain the nearly 
thermal spectra observed in most transients. If currents £11 only a tiny fraction of the 
magnetosphere, thermal photons produced in the hot surface regions have only a small 
chance of undergoing resonant scattering to populate a high-energy tail. 

Beloborodov (2009) found that the model can satisfactorily reproduce the observed 
properties of the outburst of XTE J1810-197 for a small twisted region (sin^ 9j_b ~ 0.03) 
and large twist angle (A0 ~ 1 rad). The application to other transient sources is, 
however, not without difficulties. The main problem is that the small size of the 
thermally emitting spot, and hence the limited spatial extent of the twist, can make 
the luminosity released by Ohmic dissipation too low (especially if the magnetic field 
is < 10^"^ G) to explain the observed flux, like in SGR 0501-1-4516 (Rea et ah, 2009) or 
SGR0418-1-5729, (Rea et ah, 2013a). Moreover, the relation L oc does not seem to 
be met in other sources (Rodriguez et ah, 2014). Different geometries of the j-bundle, 
not necessarily involving the polar region, may however ease the energetic requirement. 

4-2. Low-field Magnetars 

Recently, the commonly accepted picture in according to which the activity in SGRs 
and AXPs is necessarily related to a super strong-field (the ‘supercritical R’ paradigm) 
has been challenged by the discovery of two fully-fledged magnetars, SGR 0418-1-5729 
and Swift J1822.3-1606, (Rea et ah, 2010; Turolla et ah, 2011; Rea et ah, 2012a; 
Livingstone et ah, 2011b; Rea et ah, 2013a; Scholz et ah, 2012) with a dipole magnetic 
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field < 10^^ G, well within the range of ordinary radio pulsars. A third candididate, 
XMM J185246.6+003317, has been reported and awaits further monitoring (Rea et ah, 
2014a). 

SGR 0418+5729 was discovered on 2009 June 5, thanks to the detection of a couple 
of SGR-like bursts (van der Horst et ah, 2010). Follow-up observations revealed a 
previously unknown bright X-ray source at a flux level of a few times 10“^^ erg cm“^ s“^, 
pulsating at a period of 9.1 s (van der Horst et ah, 2010; Esposito et ah, 2010). Based of 
ROSAT All-Sky Survey data, the upper limit on the source flux in quiescence is of the 
order of 10“^^ erg cm“^ s“^. Despite the dense monitoring, no spin-down was detected 
during the first ~ 5 months of observations following the outburst onset. The estimated 
upper limit on the period derivative, 1.1 x 10“^^ s s“^, translates into an upper limit on 
the surface dipole magnetic held strength of 3 x 10^^ G (Esposito et ah, 2010). 

This made SGR0418+5729 the magnetar with the lowest dipole magnetic held ever 
discovered (for comparison, the previous record holder, the AXP IE 2259+586, has a 
spin-down magnetic held about twice as strong, 6 x 10^^ G). It took nearly another 
3 years of monitoring to hnally pinpoint the source spin-down rate from a coherent 
timing analysis of all the X-ray data spanning ~ 1200 days: (4 ± 1) x 10“^® s s“^, 
corresponding to R ~ 6.1 x 10^^ G and to a characteristic age Tc = P/{2P) ~ 36 Myrs 
(Esposito et ah, 2010; Rea et ah, 2010, 2013a). 

Swift J1822.3-1606 was detected on 2011 July 14, when it emitted several magnetar- 
like bursts (Livingstone et ah, 2011b; Rea et ah, 2012a, and references therein). A few 
days after the outburst onset, a new, persistent X-ray source was discovered at a hux 
level of ~ 2 X 10“^° erg cm“^ s“^, pulsating with a period of ~ 8.4 s. Gontrary to the 
case of SGR0418+5729, the source has been already detected in X-rays at a hux level 
of ~ 4 X 10“^^ erg cm“^ s“^, although its presence in two ROSAT X-ray catalogues 
passed unnoticed (Rea et ah, 2012a; Scholz et ah, 2012). Swift J1822.3-1606 was 
intensely monitored between 2011 July and 2012 August with diherent X-ray satellites 
(Livingstone et ah, 2011b; Rea et ah, 2012a; Scholz et ah, 2012). Phase coherent timing 
analyses yield spin-down rates between ~ 0.7 X 10 s s ^ and ~ 3.1 x 10 s s \ and 
a dipole magnetic held between 2.4 x 10^^ G and 5.1 x 10^^ G (Rea et ah, 2012a; Scholz et 
ah, 2012). Despite a precise measurement of P is not available as yet, any of the values 
proposed so far makes Swift J1822.3-1606 the magnetar with the second lowest dipole 
magnetic held after SGR0418+5729. Given the low value of the period derivative, the 
characteristic age of Swift J1822.3-1606 is quite long, Tc ~ 0.8 Myr, although the source 
appears not as old as SGR 0418+5729. 

The large characteristic age, the small number of detected bursts with 
comparatively low energetics and the low persistent luminosity in quiescence have been 
taken as suggestive that these are “old magnetars” approaching the end of their active 
life, in which the magnetic held has experienced substantial decay (Esposito et ah, 2010; 
Rea et ah, 2010; Turolla et ah, 2011). A key question is if, and to what extent, the 
present (internal) magnetic held is still strong enough to stress the crust and produce 
bursts/outbursts. Indeed, Turolla et al. (2011) and Rea et al. (2012a) have shown 
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Figure 19 From top left to bottom right: time evolution of P, P, Lx and Bp for 
Swift J1822.03-1606. The dashed vertical line marks the estimated age of the source; 
the gray strips in the hrst two panels show the observed values of P and P with their 
uncertainties. The model is for Btorit = 0) = 5 x 10^^ G (from Rea et al., 2012a, ©AAS. 
Reproduced with permission. A link to the original article via DOI is available in the 
electronic version ). 


that the magneto-thermal evolution (Pons et al., 2009; Vigano et ah, 2013, see §2) of 
an initially ultra-magnetized neutron star, Bp{t = 0) ~ 2 x 10^^ G, can reproduce 
the observed P, P, Bp and Lx in SGR 0418-1-5729 and Swift J1822.3-1606, for an age 
~ 1 Myr and ~ 0.5 Myr, respectively, provided that the initial internal toroidal held 
Btorit = 0) is high enough*. The evolution of the period, period derivative, dipole 
P-held and luminosity for Swift J1822.03-1606 is shown in Fig. 19. The fact that the 
characteristic age (assuming a constant B) in these two sources is largely in excess 
of that derived from magneto-thermal evolution rehects a quite general property of 
magnetar sources, for which characteristic ages are longer than those derived using 
other estimators. 

* More recent calculations of magneto-thermal evolution including the Hall drift actually show that a 
strong toroidal component develops regardless of the initial topology the the magnetic field (Vigano et 
al., 2013). 
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According to the calculations by Perna & Pons (2011), who modeled the evolution 
of the internal magnetic stresses in a magnetar, the occurrence of crustal fractures (and 
hence of bursts/outbursts) can extend to late phases (~ 10^-10® yr). Both the energetics 
and the recurrence time of the events evolve as the star ages and depend on the initial 
held. The models which successfully reproduce the properties of the two low-held sources 
imply that the two low-i? magnetars could become burst-active despite their age, with 
an expected (current) event rate of ~ 0.01-0.1 yr“^. 

Quite recently, Tiengo et al. (2013) reported the discovery of a phase-variable 
absorption feature in the X-ray spectrum of SGR 0418-1-5729. The feature is best 
detected in a 67 ks XMM observation performed on 2009 August 12, when the source 
hux was still high (5 x 10“^^ ergcm”^ s“^ in the 2-10 keV band) but is also visible in 
the RXTE and Swift data collected in the hrst two months after the outburst onset. 
The line energy is in the range ~ 1-5 keV, in XMM data, and changes sharply with 
rotational phase, by a factor of ~ 5 in one-tenth of a cycle (see Fig. 20). These seem 
to favour an interpretation in terms of a cyclotron line. If protons, e.g. contained in a 
rising hux tube close to the surface, are responsible for the line, the local value of the 
magnetic held within the baryon-loaded structure is close to 10^® G. SGR 0418-1-5729 
would then be, at the same time, the magnetar with the lowest dipole held and the 
neutron star with the largest (small-scale) held ever measured. 

4-3. Transient Radio Emission 

For a long time SGRs and AXPs were thought to be with no exceptions radio quiet, 
to the point that the lack of (pulsed) radio emission was often quoted as one of their 
dehning properties. In fact, an expanding radio nebula was detected in the aftermath 
of the Giant Flares from SGR 1806-20 (Gaensler et al., 2005b; Gameron et ah, 2005) 
and SGR 1900-1-14 (Frail, Kulkarni & Bloom, 1999), but this originated in the shocked 
material around the neutron star and not from the star magnetosphere. 

The hrst detection of pulsed radio emission from a magnetar came, rather 
unexpectedly, from the archetypal transient XTE J1810-197 (Gamilo et al., 2006), 
opening a new window in the study of magnetars. For many months, XTF J1810-197 
was the brightest radio pulsar in the Galaxy at frequencies above 20 GHz, exhibiting 
a strong variability in both the radio hux and the pulse shape on diherent timescales. 
The radio emission likely began about a year after the onset of the X-ray outburst and 
lasted a few years (Gamilo et ah, 2006, 2007a; Lazaridis et ah, 2008; Serylak et al., 2009). 
Pulsed radio emission was then discovered from the AXP IF 1547-5408 (Gamilo et al., 
2007b). This source has shown three X-ray outbursts in the past 5 years. Radio emission 
was observed in the interval between the last two events, during which it declined, to 
rise again in coincidence with the last outburst, although there was a delay of a few 
days between the onset of the X-ray outburst and the radio activity (Gamilo et al., 
2009; Burgay et ah, 2009). Another possibility is that the radio emission was blinking, 
in a way uncorrelated with the start of the X-ray outburst. PSR 1622-4950 (Levin et 
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Figure 20 The phase-dependent spectral feature in the EPIC data of SGR 0418-1-5729. 
Normalized energy versus phase image obtained by binning the EPIC source counts into 
100 phase bins and 100-eV-wide energy channels and dividing these values hrst by the 
average number of counts in the same energy bin (corresponding to the phase-averaged 
energy spectrum) and then by the relative 0.310 keV count rate in the same phase 
interval (corresponding to the pulse prohle normalized to the average count rate). The 
red line shows (for only one of the two displayed cycles) the results of a simple proton 
cyclotron model consisting of a baryon-loaded plasma loop emerging from the surface 
of a magnetar and intercepting the X-ray radiation from a small hot spot (from Tiengo 
et ah, 2013, the authors acknowledge Nature Publishing Group for reproduction). 


ah, 2010) was the hrst (and so far the only one) magnetar which was discovered thanks 
to observations in the radio band. New and archival X-ray observations have shown 
that the source is likely a transient magnetar which underwent an outburst in 2007 
(before the hrst available Chandra observation). The X-ray hux is still declining as the 
source approaches quiescence (Anderson et ah, 2012). Finally, pulsed radio emission 
has been detected from the recently discovered magnetar in the Galactic Centre, SGR 
J1745-2900 (Rea et ah, 2013b; Eatough et ah, 2013; Shannon & Johnston, 2013; Kaspi 
et ah, 2014). The source entered an outburst phase with the emission of a single 
burst on April 24 2013. The radio activity switched on 4-5 days after the outburst 
onset with similar properties to those of the other two magnetars detected at radio 
wavelengths. Despite intensive searches, pulsed radio emission was not found in other 
magnetar sources (Burgay et ah, 2006; Crawford et ah, 2007; Lazarus et ah, 2011). 

Although the sample is quite limited, a number of common features in the radio 
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emission from magnetars have started to emerge: i) association with X-ray outbursts, 
ii) the radio flux decays together with the X-ray flux but its onset is delayed, hi) marked 
variability in the radio band, and iv) flat radio spectrum. All these properties are much 
at variance with those of ordinary radio pulsars (including the high-B PSRs). Rea et 
al. (2012b) have shown that radio-loud magnetars are characterized by Lx < E in 
quiescence, much as ordinary radio pulsars, while the X-ray luminosity always exceeds 
the rotational energy loss rate in radio-silent magnetars. In fact, as indicated by Ho 
(2013), this property seems to be characteristic also of other classes of neutron stars, 
although the inverse is not true: not all sources with Lx < E emits in radio. SGR 
J1745-2900, which was not included in the original sample considered by Rea et al. 
(2012b), further conhrms this picture, having ~ 5 x 10^^ ergs“^ and Lx ^ 10^^ ergs“^ 
in quiescence (Rea et al., 2013b). This seems to point to a common origin for the 
radio emission in PSRs and in transient magnetars. In magnetars with Lx/E < 1 
particle acceleration and the subsequent ignition of the cascade process could proceed 
as in normal pulsars, and their radio emission might basically follow the same rules, 
with rotational energy driving pair creation through a cascade. The largely different 
radio properties between the two groups might result from the presence of a substantial 
toroidal component in the magnetosphere of the magnetars, contrary to the nearly 
dipolar held of PSRs. The inhuence of the large charge density required to support 
the non-potential held may also act in quenching the radio emission in the brightest 
magnetars with Lx/E > 1 (Thompson, 2008a,b), although this latter interpretation 
is hard to be reconciled with the fact that, in at least few sources (SGR J1745-2900, 
XTE J1810-197 and possibly also in IE 1547-5408 during the 2009 outburst), radio 
emission was seen while the source was bright with Lx > E. 

5. Magnetar bursts 
5.1. Burst phenomenology 

One of the hallmarks of magnetars is the emission of repeated soft gamma-ray bursts, 
and this feature played a key role in their discovery (for a nice review of the history of the 
held, see Woods & Thompson, 2006). Bursts have now been observed from 18 sources 
whose spin has also been measured (conhrming that they are indeed neutron stars) jj. 
Magnetars emit bursts in the few keV to few hundred keV energy band (hard X-ray/soft 
gamma-ray). As mentioned in the Introduction, magnetar bursts are typically grouped 
into three classes dehned primarily by energy released and their duration: short bursts 
(with energies up to 10"^^ erg), intermediate hares (energies in the range 10^^ — 10^^ erg), 
and giant hares (energies in the range 10^^ —10^® erg, see Fig. 21 for example lightcurves 
from the diherent types of burst). 

Bursting activity is highly variable. Sources can experience long periods of apparent 
quiescence (when bursting, if it occurs, is in the form of low luminosity bursts below 

tt See the Amsterdam Magnetar Burst Library, http://staff.fnwi.uva.n1/a.l.watts/magnetar/mb.html 



Magnetars: the physics behind observations 


47 






Figure 21 Examples of different types of magnetar bursts. Top left: Short bursts from 
SGR 0501+4516 recorded by Fermi/GBM (from Huppenkothen et ah, 2013). Top right: 
Burst storm from IE 1547.0-5408 recorded by Fermi/GBM, with inset showing the 
overall enhancement in emission during this event(from Kaneko et al., 2010). Lower left: 
Intermediate burst from SGR 1900+14 recorded by RXTE (from Ibrahim et ah, 2001). 
Lower right: Giant flare from SGR 1900+14 recorded by Ulysses, upper panel showing 
OTTB spectral temperature (from Hurley et ah, 1999b). The hrst three panels are 
©AAS, reproduced with permission. A link to the original articles via DOI is available 
in the electronic version. The last panel is reprinted by permission from Macmillan 
Publishers Ltd: Nature, ©1999. The authors acknowledge Nature Publishing Group 
for allowing reproduction. 

the detection threshold of our current generation of space telescopes - note also that 
there have been gaps in gamma-ray telescope coverage since the discovery of magnetars, 
due to the lack of suitable telescopes, and that sky coverage even now is never 100%), 
sporadic and highly occasional bursting, or periods of high burst activity. At their most 
dramatic, these can climax in burst storms (when several hundreds of bursts can be 
emitted over just a few hours), or rare giant flares (Fig. 21). During such outburst 
periods there may also be changes in the overall emission (luminosity, pulsed flux and 
spectrum, see Rea & Esposito 2011 for an observational review, and Pons & Rea 2012 
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for a recent theoretical study) and the timing behaviour (spin-down rate and glitching, 
see for example Woods et ah 2001, 2002, 2003, 2007; Dib & Kaspi 2014). There is 
no clear dependence of short/intermediate bursting activity on the dipolar magnetic 
held strength of the source. The most energetic giant hares come from three sources 
whose inferred dipole helds are amongst the highest (a few times 10^^ G up to ~ 10^^ 
G). However bursts have also been seen from sources with helds apparently below the 
quantum critical limit, such as Swift J1822.3-1606 with a held ~ 1.35 x 10^^ G (Scholz, 
Kaspi & Gumming, 2014) and SGR 0418-1-5729 with held 6 x 10^^ G (Rea et al., 2010)//. 

5.1.1. Short bursts Short bursts are the most common, with huences in the range 
10“^^ — 10“^ erg/cm^/s, which for the assumed distances implies isotropic energies in 
the range 10^® — 10"^^ erg (Gogii§ et al., 1999, 2000; Woods et al., 1999b; Gavriil, Kaspi 
& Woods, 2004; Kumar, Ibrahim & Sah-Harb, 2010; Scholz & Kaspi, 2011; Lin et al., 
2013). Peak luminosities are in the range 10^® — 10"^^ erg/s, extending to well above the 
standard Eddington luminosity of ~ 2 x 10^® erg/s for a non-magnetic neutron star. The 
lowest luminosity bursts recorded are at the sensitivity limit of our current generation 
of detectors. For sources that have shown sufficient numbers of bursts for meaningful 
statistical analysis (SGR 1806-20, SGR 1900+14, SGR 1627-41, IE 2259+286, SGR 
0501+5416 and IE 1547.0-5408), burst huences are distributed as a power law, dN = 
E”''' dE with 7 ~ 1.4 — 2.0 (Gheng et ah, 1996; Aptekar et ah, 2001; Woods et al., 
1999b; Gogii§ et ah, 1999, 2000, 2001; Gavriil, Kaspi & Woods, 2004; Kumar, Ibrahim 
& Sah-Harb, 2010; Savchenko et ah, 2010; Lin et al., 2011a; Scholz & Kaspi, 2011; 
Prieskorn & Kaaret, 2012; van der Horst et al., 2012). Burst durations are ~ 0.01 — Is 
and are distributed lognormally with a peak at ~ 0.1s, less than the rotational period 
of the star. Duration is correlated with huence. Bursts from the sources that have 
shown insufficient events for full statistical analysis are consistent with this picture. 
Lightcurves are extremely variable in shape: the rise is in general faster than the decay 
(Gogii§ et ah, 2001; Gavriil, Kaspi & Woods, 2004), but bursts can be multi-peaked, and 
no simple phenomenological model has yet been found that would ht the morphologies of 
the diherent burst lightcurves. Some short bursts (classihed according to their fluence) 
from hve sources have extended faint tails of ~ 100 — 1000 s in duration, leading to 
an overall energy release that can exceed that of the original spike (Woods et ah, 2005; 
Gavriil, Kaspi & Woods, 2004, 2002, 2006; Dib, Kaspi & Gavriil, 2009; An et ah, 2014a; 
Gavriil, Dib & Kaspi, 2011; Mereghetti et ah, 2009; Savchenko et ah, 2010; Scholz & 
Kaspi, 2011). The tails appear to be pulsed at the rotational frequency. 

From the sources with large samples, wait times (for bursts above the detection 
threshold, in periods of continuous telescope coverage/) form a lognormal distribution 

// Although note that the dipole field strength as estimated from spin down provides only a lower limit 
on the magnetic field strength, and there are no constraints on the strength of either higher order 
poloidal components or toroidal components. 

/ The recent detection of short bursts from IE 1547.0-5408, in the VLF radio band due to the 
ionospheric disturbance that occurs as the incident gamma-rays ionize the Earth’s atmosphere opens 
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with a peak at ~ 100s. During the active period of SGR 1900+14 in 1998, for example, 
Gogu§ et ah (1999) found wait times ranging from less than Is to more than lOOOs:};, 
and far longer wait times are clearly possible: SGR 1900+14, for example, has also had 
quiescent periods ~ years in duration. The shortest wait times observed are comparable 
to the durations of individual bursts, such that the distinction between single bursts 
and multi-peaked events is not clear. There appears to be no correlation between 
burst intensity and wait time to the following burst (Laros et ah, 1987; Gogu§ et ah, 
2000; Gavriil, Kaspi & Woods, 2004; Savchenko et ah, 2010). Of particular note in the 
discussion of wait times are burst storms, periods of unusually high short burst activity 
(which may include some intermediate flares), in which tens to hundreds of bursts occur 
over only a few hours on top of an overall rise in emission that can be strongly pulsed 
at the rotational phase (see for example Hurley et ah, 1999a; Israel et ah, 2008; Gavriil, 
Kaspi & Woods, 2004; Mereghetti et ah, 2009; Kaneko et ah, 2010; Savchenko et ah, 
2010, for the cases SGR 1900+14, IE 2259+286, and IE 1547.0-5408). There have also 
been efforts to determine whether the occurrence of bursts correlates with rotational 
phase: here the evidence is mixed. For IE 1048.1-5937, IE 2259+286 and XTE J1810- 
197 bursts do seem to occur preferentially at rotational pulse maxima (Gavriil, Kaspi 
& Woods, 2002, 2004; Woods et ah, 2005). However for SGR 1806-20, SGR 1900+14, 
SGR 1627-41, 4U 0142+61 and IE 1547.0-5408 no such correlation is found (Palmer, 
1999, 2002; Woods et ah, 1999b; Gavriil, Dib & Kaspi, 2011; Savchenko et ah, 2010; 
Scholz & Kaspi, 2011; Lin et ah, 2012a). 

5.1.2. Intermediate flares Intermediate flares, with (isotropic) energies in the range 
~ 10"^^ — 10^^ erg, and peak luminosities that exceed the non-magnetic Eddington limit, 
have been seen from SGR 1627-41, SGR 1900+14, SGR 1806-20, and IE 1547.0-5408. 
The primary bursts appear to be brighter and slightly longer (durations ~ 0.5s up to 
a few s) versions of the short bursts. Morphologies, however, are varied. Some have a 
clear decay and an abrupt end (Mazets et ah, 1999a,b; Olive et ah, 2004; Israel et ah, 
2008). In others the initial burst is followed by a extended decaying tail that can last for 
up to several thousand seconds, but contains less than ~ 2 % of the energy released in 
the initial peak (Ibrahim et ah, 2001; Lenters et ah, 2003; Esposito et ah, 2007b; Gogu§ 
et ah, 2011b). The tails are pulsed at the rotational period of the star: in some cases 
the pulsed amplitude rises dramatically (Ibrahim et ah, 2001; Lenters et ah, 2003), as 
seen during some burst storms (Kaneko et ah, 2010); whilst in others no change is seen 
(Gogu§ et ah, 2011b). The pulsations appear for the most part to be phase-aligned with 
the pre-burst pulsations, but there are some occasional exceptions (Guidorzi et ah, 2004; 
Gogu§ et ah, 2011b). There has also been one burst with a decaying pulsed tail, where 
the sharp initial burst peak appears to be absent, leading to a rather slow rise time 
~ 10s (Kouveliotou et ah, 2001; Guidorzi et ah, 2004). Intermediate flares sometimes 

up the possibility of using the VLF band to obtain a more complete coverage of the waiting time 
distribution with being dependent on sky coverage of space telescopes (Tanaka et ah, 2010). 
f The longest wait times in this study were set by the length of the observing window. 
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have short precursors (Ibrahim et ah, 2001; Gogii§ et ah, 2011b), and during some events 
short bursts are seen during the extended tail (Ibrahim et ah, 2001; Lenters et al., 2003; 
Gogu§ et ah, 2011b). 

5.1.3. Giant flares The most energetic bursts, the giant flares, are extremely rare. Only 
three have ever been seen, in 1979, 1998 and 2004, each from a different magnetar (SGR 
0526-66, SGR 1900-1-14, and SGR 1806-20). The total energy released, if the emission 
is isotropic and assuming reliable estimates of distance, is in the range 10*^^ — 10^^ erg 
(Fenimore et ah, 1996; Feroci et ah, 2001; Palmer et al., 2005). The overall properties 
of the three giant flares are, despite the differences in energy, very similar (in marked 
contrast with the heterogeneity of the intermediate flares). They have a very bright 
initial peak, followed by an extended decaying tail with a duration of several hundred 
seconds that is strongly pulsed at the rotational frequency of the starf (see e.g. Mazets 
et al., 1979a; Hurley et al., 1999b; Hurley et ah, 2005). 

The initial peaks, which can have rise times as short as ~ 1 ms, last ~ 0.1 — 1 
s and are very hard. Luminosities reach up to 10^'^ erg/s (Hurley et ah, 2005), which 
causes substantial dead time and pile up effects in space telescopes. This renders reliable 
spectral modelling very difficult, however the spectrum is very hard, with emission being 
detected up to 2 MeV (Mazets et ah, 1979a; Hurley et ah, 1999b). The initial peaks are 
strongly variable, on timescales as short as a few ms (Barat et al., 1979; Hurley et al., 
1999b; Terasawa et al., 2005; Schwartz et al., 2005). Both the SGR 1900-1-14 and SGR 
1806-20 giant flares were observed to have precursorsf. For SGR 1900-1-14 the precursor 
resembled a normal short burst, and occurred < 1 s before the giant flare (Mazets et 
al., 1999c). For SGR 1806-20 the precursor was flat-topped, with an energy release that 
puts it in the intermediate flare class, and occurred 142 s prior to the giant flare (Hurley 
et ah, 2005). A discussion of whether the apparent precursors are in fact genuinely 
causally connected to the giant flares was presented by Gill & Heyl (2010). 

The energies emitted in the tails of the three giant flares have been similar (~ 10^"^ 
erg). This is 1 — 2% of the energy released in the initial peak of the 2004 giant flare: for 
the two earlier giant flares the energies released in initial peak and tail were comparable. 
The overall envelope of the tails (averaged over rotational phase) decays smoothly as a 
power law, coming to an abrupt end after several hundred seconds. The pulsations in 
the SGR 0526-66 giant flare tail were seen immediately after the initial peak; for the 
other two giant flares they appeared only a few tens of seconds later. Pulse prohles 
can evolve during the tails, in the case of SGR 1900-1-14 simplifying quite dramatically 
(Feroci et ah, 2001), with evolution being much more minimal for SGR 1806-20 (Palmer 

f Note that none of the giant flares were caught during pointed observations: they are so rare but so 
bright that most are seen off-axis. This means that the sensitivity to late time weak emission is much 
less than for some of the intermediate flares, which were observed during pointed observations. One 
should bear this in mind when comparing the apparent durations. 

t A precursor with the properties observed for SGR 1900-1-14 giant flare would not have been detectable 
for the SGR 0526-66 giant flare given the instrumentation at the time (Gill & Heyl, 2010). 
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et al., 2005; Mereghetti et al., 2005c; Boggs et al., 2007; Xing & Yu, 2011). The giant 
flares are so strong that they have a detectable effect on the Earth’s electromagnetic 
held (Mandea & Balasis, 2006) and ionosphere (Inan et ah, 1999, 2007; Tanaka et ah, 
2008), with even the rotational pulsations being clearly visible, and this has been used 
to put lower limits on the strength of low energy (< 10 keV) emission from the burst 
(unaffected by satellite dead time issues). 

Radio afterglows were detected after both the SGR 1900+14 and SGR 1806-20 
giant flares (Frail, Kulkarni & Bloom, 1999; Gameron et ah, 2005; Gaensler et ah, 2005b; 
Gelfand et ah, 2005; Taylor et ah, 2005; Fender et ah, 2006). The amount of energy in 
the radio afterglow is much less than that emitted in the gamma-rays. This is different 
from what is observed in gamma ray bursts, where the ratio between the two types 
of emission is of order unity, such that the lower energy and longer duration emission 
is associated with re-processing of the gamma-ray energy by the surrounding material. 
The radio afterglow from the magnetar giant flares is linearly polarized, implying that 
it is caused by electron synchrotron emission, and is observed to expand over time. Its 
generation requires an ejection of relativistic particles and magnetic fields (in the form 
of a “plasmon”, which expands and cools) by the burst process (see Section 5.3). 

5.2. Burst trigger mechanisms 

Rapid magnetic field reconfiguration is assumed to be an integral part of the bursts: as 
we will discuss in Sec. 5.3, the gamma-ray emission is assumed to come from particles 
accelerated by rapid held change. Thus slow magnetic evolution builds up stresses in 
the system, some of which are released catastrophically in bursts, which must either be 
driven by or result in rapid magnetic field reconfiguration. However the precise trigger 
mechanism, and the role of the magnetic held within it (if any), remains unclear. Three 
main locations (and associated families of instabilities) have been considered for the 
trigger mechanism: below we review each in turn. 

The first option, as suggested by Thompson & Duncan (1995), is that the magnetic 
held evolves into an unstable configuration within the liquid core of the star (Markey, 
1973; Wright, 1973; Tayler, 1973; Flowers & Ruderman, 1977), which is then susceptible 
to a large-scale magnetohydrodynamical instability (Lander & Jones, 2011a; Giolfl et ah, 
2011; Kiuchi, Yoshida & Shibata, 2011; Giolfl & Rezzolla, 2012). This would develop on 
the Alfven crossing time of the core, which is ~ 0.1s for a 10^^ G interior field (r ~ R/va, 
where the Alfven speed va is given in Equation 16 later in this paper) and hence broadly 
compatible with the durations of both the normal bursts and the initial peaks of the 
giant flares. 

Whether such an unstable state could develop is open to question, since if the core 
is superconducting, the high conductivity should facilitate swift reconfiguration and 
hence prevent the formation of such an unstable state. However the configuration could 
be stabilized by currents in the crust, or superconductivity may be suppressed if the 
core held exceeds ~ 10^®G. Such rapid internal magnetic reconfiguration would inject 
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an Alfven pulse into the magnetosphere that would then generate the observed burst 
emission. The energy available to power the burst in this scenario, since it originates in 
the core, is more than sufficient to power even a giant flare. Some internal heat release 
may also be expected due to dissipation in core and crust as the instability proceeds. 

The second option is that the decay of the core held places magnetic stresses on 
the solid crust of the star (an ionic lattice to which the held in the crust is locked). 
The crust can deform elastically to accommodate this up to a certain point, then 
ruptures catastrophically once its breaking strain is exceeded (Thompson & Duncan, 
1995; Thompson & Duncan, 2001). The stored energy released in this scenario would 
come from the crust and possibly also the core: although initially it was thought that the 
crust alone could not store enough elastic energy to power the giant flares (Thompson 
& Duncan, 1995), the latest molecular dynamics simulations predict a higher breaking 
strain, indicating that this may be feasible after all (Horowitz & Kadau, 2009; Hoffman 
& Heyl, 2012). Studies are now underway that aim to determine how and where stresses 
would build up in the crust as a result of core held evolution, the goal being to determine 
how often, and where (location including depth) the crust is most likely to fail (Perna 
& Pons, 2011; Pons & Perna, 2011; Beloborodov & Levin, 2014; Lander et al., 2015). 
Recent calculations by Lander et al. (2015) of the strain induced in a crust by a changing 
magnetic held configuration And a characteristic burst energy 


_^_~ 0.25 ^ ^ 
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where (Jmax is the breaking strain of the crust (which could be as high as 0.1; Horowitz 
& Kadau, 2009), d is the depth at which the crust ruptures, Rc the crust thickness, R is 
the star radius, and I the rupture length. The characteristic local held strength related 
to crust breaking (from the same study) is given by 

Bw=.k = 2.4xlO“G(^)‘''V (12) 

When the crust does rupture it must do so by rapid plastic deformation, not via 
brittle fracture, due to the impossibility of opening up voids under the conditions of 
extreme pressure that pervade in neutron star crusts (Jones, 2003). The role of the 
magnetic held during crust rupture however is not clear: Levin & Lyutikov (2012) have 
argued that under some circumstances the held deformation induced by an incipient 
rupture may act as a brake on its propagation. Shear wave timescales (which control 
the timescale on which the crust ruptures, r ~ nR/vg, where the shear speed Vg is 
given by Equation 15 later in this paper) are compatible with those observed in bursts 
and flares (see for example Schwartz et al., 2005), injecting an Alfven pulse into the 
magnetosphere. However the transfer of energy into the external magnetosphere may 
be slowed by a large impedance mismatch at the crust-magnetosphere boundary (Link, 
2014). Crustal rupture would most likely lead to local heating as well. 

The final possibility is that the core and crust evolve smoothly, and that stress builds 
up instead in the magnetosphere. Stress release is then envisaged as taking place via 
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a plasma instability involving spontaneous magnetic reconnection (see Uzdensky, 2011, 
for a review). A number of studies have looked at how the external magnetosphere 
might respond to the expulsion of magnetic helicity due to the decaying core field, and 
have found that the development of unstable conhgurations with strong magnetic shear, 
that might be prone to reconnection instabilities, is feasible (Thompson, Lyutikov & 
Kulkarni, 2002; Beloborodov, 2009; Parfrey, Beloborodov & Hui, 2012; Parfrey et ah, 
2013). The resulting explosive reconnection event would progress on the Alfven crossing 
time in the magnetosphere, which is < 0.01 s (since in the magnetosphere va ~ c), 
and the energy that can be released is more than sufficient to power even the giant 
flares. Specihc instabilities that have been considered primarily for the giant flares 
(driven by earlier concerns about the ability of crust ruptures to release enough energy) 
are the relativistic tearing mode (Lyutikov, 2003, 2006; Komissarov et ah, 2007) and 
collisionless Hall reconnection mediated by emission from precursor bursts (Gill & Heyl, 
2010), the precursors presumably being triggered in this scenario by another mechanism. 
Similarities between the giant flares and reconnection driven coronal mass and flux tube 
ejection events in solar physics have also been explored in some depth (Masada et ah, 
2010; Yu, 2011, 2012; Yu & Huang, 2013; Huang & Yu, 2014a,b; Meng et ah, 2014). 
Instabilities may also arise from the interaction of MHD waves with the vacuum in helds 
above the quantum critical limit (where QED effects are important, Heyl & Hernquist, 
2005). All of these mechanisms would lead directly to particle acceleration and radiation 
in the magnetosphere, with possible crustal heating via particles impacting the surface. 

At present it is by no means clear which of the various mechanisms are in operation, 
and given the diversity of burst properties (Sec. 5.1) more than one may be in operation. 
Timescales, as discussed above, appear to be roughly compatible with either burst rise 
times or durations (assuming that the emission process timescales reflect those of the 
trigger mechanism). Moreover both starquakes and magnetospheric reconnection could 
in principle explain the power law distribution of fluences (which is often taken as 
evidence for Self-Organised Criticality, see for example Aschwanden et ah, 2014). Serious 
efforts are now being made to simulate the build up of stress and the development of 
instabilities, but in order to allow meaningful tests of the data, consideration must be 
given to how these various triggers connect to the emission that we see. We discuss this 
in the next section (Sec. 5.3). 

5.3. Burst emission processes 

5.3.1. Sources of emission Magnetar bursts are complex, with varied spectra and 
morphologies. The emission process for all bursts is generally assumed to be started by 
rapid rearrangement of the magnetic held (resulting from one of the trigger mechanisms 
discussed in Sec. 5.2), possibly involving either induced or spontaneous magnetic 
reconnection. This accelerates charged particles with ensuing gamma-ray emission, 
since the rapid acceleration of electrons in a strong curved held leads to a cascade of pair 
creation and gamma-rays (Sturrock, Harding & Daugherty, 1989). To obtain the hardest 
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emission, there must be very low contamination by baryonic material (since scattering 
would lead to softening). Fully self-consistent models of the emission resulting from the 
various proposed trigger mechanisms do not yet exist. Despite this, the rise timescale 
of the gamma-ray emission has frequently been used as a key piece of evidence to argue 
for a particular trigger mechanism, as discussed in the previous section. However it is 
not clear that this is warranted: details of the gamma-ray emission process may in fact 
completely obscure the timescales associated with the original trigger mechanism (see 
for example Hoshino & Lyubarsky, 2012). 

To obtain the radio afterglow seen in the giant flares, it is necessary to postulate 
the ejection of a plasmoid of magnetic fields and trapped shocked plasma, that gradually 
cools (Sec. 5.1). Such plasmoid ejection is a natural and expected consequence of a large- 
scale reconnection event in the magnetosphere (see Sec. 5.2 and for example Lyutikov, 
2006). 

The initial spike of magnetar flares may also lead to radio emission (explored for 
example in Lyutikov, 2002). More recently it has been suggested that the interaction 
between strongly magnetized relativistic ejecta (expelled by the initial spike of giant 
flares) and the surrounding wind nebula might be responsible for extragalactic Fast 
Radio Bursts (Popov & Postnov, 2007, 2013; Thornton et ah, 2013; Lyubarsky, 2014). 
At present this remains speculative. 

If the local energy generation rate is high enough, as explained above, it will lead to 
copious production of electron-positron pairs and gamma-rays. If this occurs in a closed 
held line region, where the charged pairs cannot cross magnetic held lines, they become 
trapped. As density increases, so does optical thickness, trapping the photons as well 
and leading to rapid thermalization (Thompson & Duncan, 1995). The held necessary 
to conhne the plasma can be estimated by requiring that magnetic pressure exceed the 
pressure of the radiation and the pairs at the outer boundary of the hreball, yielding 
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(Thompson & Duncan, 2001), where AR is the characteristic size of the hreball and R 
the neutron star radius. 

Such a trapped pair plasma hreball would then cool and contract due to radiative 
dihusion from a thin surface layer, with the bulk of the radiation leakage occurring close 
to the stellar surface since this is where the held is strongest and scattering the most 
suppressed. The opacity at the surface will be dominated by the electron-ion plasma 
ablated from the neutron star surface (especially if the emergent hux is close to the 
magnetic Eddington limit) which form the photosphere. This trapping would prolong 
the emission from the burst by acting as a reservoir for the energy. Fireball formation is 
not exclusively linked to any one trigger mechanism (Thompson & Duncan, 1995; Heyl 
& Hernquist, 2005). 

The hreball model has been very successful at explaining the later decaying tail 
phase of giant hares (Thompson & Duncan, 2001). A cooling hreball trapped on 
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closed field lines should have a luminosity L whose time-dependence is described by 
the following function 


L{t) = L(0) 
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where the cooling luminosity is assumed to vary as a power of the remaining hreball 
energy, L oc (Thompson & Duncan, 2001). This proves to be a good £t to giant flare 
tail data, with Tevap of order a few hundred seconds, and a value of a that is close to that 
expected for a spherical fireball of uniform temperature (Feroci et al., 2001). Spectral 
fitting indicates that the emitting area falls while the temperature of the radiation 
remains roughly constant at the level expected for the photosphere of a trapped fireball 
in a magnetic field in excess of the quantum critical field (Thompson & Duncan, 1995; 
Feroci et al., 2001). The observed photospheric temperature (~ 20 — 30 keV) is lower 
than the inferred temperature in the core of the hreball, which is ~ 100 keV (Thompson 
& Duncan, 1995; Thompson & Duncan, 2001). Thompson & Duncan (1995) argue 
that this is an intrinsic property of the way various processes act to preserve thermal 
equilibrium in the fireball, and the way that radiation gradually escapes. However 
other processes such as photon splitting as the emitted radiation propagates through 
the strong magnetic held may also be important (Baring, 1995). The beaming of 
radiation as it leaks from the base of the fireball and streams along held lines provides a 
simple explanation for the strong rotational pulses seen in the giant hare tails (Fig. 22). 
Whether fireballs form in the smaller bursts is still not clear (energy release may not 
occur at a fast enough rate (Gogii§ et al., 2001), although the similarity of the spectra 
of the short bursts to the spectra in the tails of the giant hares suggests that there may 
be a link. 

Some of the energy released during the burst is likely to excite vibrations, either of 
the star (crust/core, see Sec. 5.4) or in the form of Alfven waves in the magnetosphere. 
This too can act as a store for energy that is then radiated on longer timescales. If the 
vibration rate is slow enough, this can act as a source of ongoing excitation that forms 
an extended pair corona (that obscures and scatters the radiation from the trapped 
hreball) from which the radiation emerges isotropically (Feroci et ah, 2001; Thompson 
& Duncan, 2001). The presence of such an extended pair corona has been invoked to 
explain the smooth emission immediately after the peak of the giant hares, which then 
clears over 30-40 s to reveal the strongly beamed rotational pulse emanating from the 
trapped hreball beneath (Feroci et ah, 2001). 

A burst may also have a thermal component of emission that is produced by residual 
heat of crust rupturing (Lyubarsky, Eichler & Thompson, 2002; Kouveliotou et al., 
2003), extreme heating and possible melting of the crust immediately underneath a 
trapped hreball (Thompson & Duncan, 1995), or bombardment of the stellar surface 
by magnetospheric particles (Lyutikov, 2006; Beloborodov, 2009). Such thermal 
components, particularly deep crustal heating, are one possible explanation for the 
presence of both the additional pulsed components seen after some bursts (short and 
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Figure 22 Left: Figure from Thompson & Duncan (2001), to illustrate how radiation 
escapes from a trapped photon-pair plasma fireball. Scattering opacities are strongly 
polarization dependent, with most radiation escaping from the hreball in the E-mode. 
This radiation is then collimated along open magnetic held lines (due to the magnetic 
held dependence of the scattering opacties), forming a beam that is then modulated 
by the star’s rotation to give rise to a rotational pulse (from Thompson & Duncan, 
2001, ©AAS. Reproduced with permission. A link to the original article via DOI is 
available in the electronic version). Right; Figure from Lyubarsky (2002) showing how 
the photospheric spectrum from a hreball with temperature 15keV (dashed line) would 
be modihed by various high held radiation processes: cyclotron absorption, for a 6 x 10^"^ 
G held, and photon splitting (from Lyubarsky, 2002, with OUP permission). 


intermediate) - as a localised hotspot - and the must longer decaying afterglows seen 
after burst active periods, intermediate hares and giant hares (Kouveliotou et ah, 2003; 
Feroci et ah, 2003). 

The location of the emitting regions is further complicated by the fact that the 
luminosity of the bright bursts may exceed the relevant Eddington limit, leading to 
photospheric expansion and ejection of material as radiation pressure overwhelms the 
gravitational force (see for example Thompson & Duncan, 1995; Watts et ah, 2010). 
Magnetar bursts can easily exceed the non-magnetic Eddington limit. However strong 
magnetic helds suppress scattering opacities, increasing the Eddington luminosity even 
before magnetic conhnement ehects - which can increase the limit still further - are taken 
into account (Paczynski, 1992; Thompson & Duncan, 1995; Miller, 1995; van Putten et 
ah, 2013). 

5.3.2. Radiative transfer processes Short burst spectra, in an era of improved 
broadband coverage, are typically well ht as either two blackbodies (2BB, with 
temperatures ~ 5 keV and ~ 15 keV) or using a Comptonization (power law with a 
high energy cutoff) model (Feroci et ah, 2004; Olive et ah, 2004; Nakagawa et ah, 2007; 
Israel et ah, 2008; Esposito et ah, 2008; Lin et ah, 2011a; Scholz & Kaspi, 2011; van 
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der Horst et al., 2012; Lin et al., 2012a). There is a sharp correlation between radius 
and temperature of the blackbodies in the 2BB hts, with the softer BB component 
saturating before the harder one as burst fluence increases. For some sources bursts 
tend to harden with increasing fluence, whereas for others they soften (Gogii§ et ah, 
2001; Gavriil, Kaspi & Woods, 2004; Gotz et ah, 2004; Kumar, Ibrahim & Safi-Harb, 
2010; Savchenko et ah, 2010; Scholz & Kaspi, 2011; van der Horst et ah, 2012). Earlier 
burst papers tend to use an Optically Thin Thermal Bremsstrahlung (OTTB) model 
(Gogii§ et ah, 1999; Woods et ah, 1999b; Gogii§ et ah, 2000; Aptekar et ah, 2001), 
which fit the data well above 15 keV although they overpredict the flux of photons at 
low energies (Fenimore, Laros & Ulmer, 1994; Feroci et ah, 2004). Despite this, OTTB 
hts are often included in more recent analysis, to allow comparison with earlier studies. 
OTTB temperatures are typically in the range 20-40 keV, and in general no emission 
is seen from short bursts above 150-200 keV. There have however been a handful of 
events from SGR 1900-1-14 with a much harder spectrum, and emission extending up to 
500 keV (Woods et ah, 1999d). There have also now been studies exploring the softest 
part of the burst spectrum, using data from XMM-Newton, where the burst spectra 
appear to be well ht with a more physically-motivated model comprising a modihed 
blackbody plus resonant cyclotron scattering Lin et al. (2012b, 2013, and see Sec. 5.3). 
There have also been strong efforts to search for spectral lines in magnetar bursts. For 
a long time the only reported detections, from 5 keV to 13-14 keV, in bursts from 
SGR 1806-20, XTE J1810-197, 4U 0142-1-61, and 1E1048.1-5937 came from RXTE data 
(Ibrahim et ah, 2002; Woods et ah, 2005; Gavriil, Kaspi & Woods, 2002, 2006; Gavriil, 
Dib & Kaspi, 2011). However recently a similar feature has been detected in bursts from 
1E1048.1-5937 observed by NuSTAR (An et ah, 2014a), increasing conhdence that they 
are indeed intrinsic to the bursts. The line energy is close to that expected for the proton 
cyclotron line given the inferred magnetic held strength. In addition to time-integrated 
spectra, data quality are now sufficiently good that it is possible to do time-resolved 
spectroscopy. These studies indicate that although the best ht spectral model remains 
the same during individual bursts, the parameters can evolve (Israel et ah, 2008; Lin et 
ah, 2011a; Younes et ah, 2014). However between bursting episodes, the best ht model 
for individual sources may change (von Kienlin et ah, 2012). 

The spectra of the initial spikes of intermediate hares are similar to the short bursts 
(Mazets et ah, 1999a,b; Olive et ah, 2004; Israel et ah, 2008). The spectrum of the 
extended decaying tails is however diherent (in contrast with the tails seen after some 
short bursts, where the peak and tail have similar spectra). Tail spectra for intermediate 
hares are well ht by a BB, possibly with an additional power law component, and the 
emission softens during the tail (Ibrahim et ah, 2001; Tenters et ah, 2003; Esposito et ah, 
2007b; Gogii§ et ah, 2011b). 

Eor the giant hares, reliable spectral modelling in the initial spike is complicated 
enormously by dead time and pile up (Fenimore et ah, 1996; Mazets et ah, 1999c). 
However OTTB models or quasi-BB models yield spectral temperatures in the range 
200-300 keV, and emission has been detected up to 2 MeV (Mazets et ah, 1979a; Hurley 
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et al., 1999b; Hurley et al., 2005). The spectrum of the emission in the tail is very 
similar to that of the short bursts, with OTTB temperatures ~ 10 — 30 keV (Fenimore 
et ah, 1981; Hurley et ah, 1999b) that vary with rotational phase. Other spectral models 
such as BB or 2BB, possibly with a power law, also provide a good £t to the data. A 
significant hard (> 1 MeV) component was seen during the tail and subsequent afterglow 
from the SGR 1806-20 giant flare (Mereghetti et ah, 2005c; Frederiks et ah, 2007; Boggs 
et ah, 2007). 

Although the spectral models described above provide a reasonable fit for the data, 
they are not based on physical models that take into account all of the scattering and 
resonant processes known to be important in such strong magnetic fields. Any thermal 
emission, for example, as might be expected from a trapped fireball (Thompson & 
Duncan, 1995; Thompson & Duncan, 2001), would be strongly modified, see Fig. 22. 
Lower energy photons scatter less and can hence escape from deeper, hotter parts of the 
atmosphere. The radiation at low energies should thus exceed that expected for simple 
blackbody emission (Ulmer, 1994; Lyubarsky, 2002). Photon splitting and merging will 
also be important in modifying the spectrum (Miller, 1995; Baring, 1995; Thompson & 
Duncan, 2001) at energies above around 30 keV, and resonant cyclotron scattering (RCS, 
see Sec. 3.3.1) will also be important. Efforts to fit burst spectra using more physical 
models, or to intepret the phenomenological models in terms of physical parameters, 
are however very rare. Israel et al. (2008) (also Kumar, Ibrahim & Safi-Harb (2010)) 
suggested that the two blackbodies in the 2BB model fits might be the photospheres 
associated with the different poiarization modes, although theoretical calculations of 
the properties of the two photospheres do not match those inferred from the fits (van 
Putten et ah, 2013). Lin et al. (2011a) attempted to interpret the parameters of the 
Comptonization and 2BB model fits in terms of a population of coronal electrons 
scattering surface emission (for example from a fireball), see also van der Horst et 
al. (2012) and Yonnes et al. (2014). More recently, Lin et al. (2012b, 2013) made 
an effort to fit soft burst emission using the modified biackbody model developed by 
Lyubarsky (2002), augmented to include effects of RCS. The fact that physical model 
interpretations are still so scarce, however, emphasizes the hnge nncertainty in terms of 
the location of the emission mechanisms, how they form, and how the released energy 
is partitioned between them. 

5.4- Burst seismology 

Asteroseismology is a precision technique for the study of stellar interiors, and it is 
magnetars that have opened up this field for neutron stars. This began when Qnasi- 
Periodic Oscillations (QPOs) in the hard X-ray emission were fonnd in the tails of the 
giant flares from the magnetars SGR 1806-20 (Israel et ah, 2005; Watts & Strohmayer, 
2006; Strohmayer & Watts, 2006) and SGR 1900-1-14 (Strohmayer & Watts, 2005). In 
the tail of the SGR 1806-20 giant flare (Fig. 23) there were several QPOs in the range 
18-150 Hz, and two isolated higher frequency signals at 625 Hz and 1840 Hz. The QPOs 
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Figure 23 Figure from Strohmayer & Watts (2006) showing the two strongest QPOs 
detected in the tail of the SGR 1806-20 giant flare. Top left: RXTE lightcurve of 
the event. Lower left: Rotational pulse during this time: the power spectra shown are 
computed using the segments enclosed by the dashed lines. Right: Power spectrum made 
by averaging nine 3 s segments from the time interval marked by dashed lines in the top 
left panel. The 92 Hz and 625 Hz QPOs are clearly visible, and the inset illustrates the 
signihcance of the 625 Hz feature (from Strohmayer & Watts, 2006, ©AAS. Reproduced 
with permission. A link to the original article via DOI is available in the electronic 
version) 


detected in the tail of the giant flare from SGR 1900-1-14 had frequencies in the range 
28-155 Hz. Widths (FWHM) were in the range 1-20 Hz, with fractional amplitudes up 
to ~ 20 % rms that are strongly rotational phase-dependent. 

The idea that giant flares might excite global seismic vibrations was hrst predicted 
by Duncan (1998), and this is the most plausible explanation that has yet been advanced 
to explain the QPOs (Israel et ah, 2005). If this interpretation is correct, such vibrations 
offer an unprecedented opportunity to constrain the interior held strength and geometry 
(something that is very hard to measure directly), and also perhaps the dense matter 
equation of state (Samuelsson & Andersson, 2007; Watts & Reddy, 2007). In order to 
do this, however, the modes must be correctly identihed. 

The QPOs were initially tentatively identihed with torsional shear modes of the 
neutron star crust and torsional Alfven modes of the highly magnetized huid core. 
These identihcations were based on the expected mode frequencies, which are set by 
both the size of the resonant volume and the relevant wave speed. For crustal shear 
modes, the appropriate speed is the shear speed Vg = {p,s /where ps is the shear 
modulus and p the density. The shear modulus is of the order of the Goulomb potential 
energy ~ Z'^e^/r per unit volume r^, where r ~ [p /is the inter-ion spacing, 
while Z and A are the effective atomic number and mass number, respectively, of the 
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ions in the crust. Using the shear modulus computed by Strohmayer et al. (1991) and 
scaling by typical values for the inner crust (Douchin & Haensel, 2001), the shear velocity 
as shown by Piro (2005) is: 
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where is the fraction of neutrons. This yields a rough estimate for the frequency 
for the fundamental crustal shear mode of z/ ~ Vg/^nR = 18 (10 km/i?) Hz. Full mode 
calculations hud similar values, but with additional dependencies on the mass and radius 
of the star due to relativistic effects (see for example Samuelsson & Andersson 2007), 
and it is this dependence that makes the modes potentially powerful diagnostics of the 
dense matter equation of state (Lattimer & Prakash, 2007). Many of the lower QPO 
frequencies could be explained as angular harmonics with no radial nodes, whilst the 
two highest frequencies in the SGR 1806-20 giant flare were identihed as radial overtones 
of these crustal modes. 

For torsional Alfven modes of the core, the appropriate wave speed is the Alfven 
speed va = B /where B is the magnetic held strength, giving 
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This yields a very rough estimate for the frequency of the fundamental torsional Alfven 
mode of z/ ~ va/'^R = 25 (10 km/R) Hz (Thompson & Duncan, 2001). Note however 
that the value of the held strength B in magnetar cores is highly uncertain, as is the 
appropriate value of the density p. In principle only the charged component (~ 5- 
10% of the core mass) should participate in Alfven oscillations, reducing p, however 
there are mechanisms associated with superhuidity and superconductivity that can 
couple the charged and neutral components, leading to additional mass-loading. As 
above, full mode calculations that take into account relativistic ehects lead to additional 
dependencies on neutron star mass and radius (see for example Sotani et al. 2008). It 
should also be noted that the Alfven modes constitute continua rather than a set of 
discrete frequencies, since the held lines within the core have a continuum of lengths. The 
observed QPOs would then be associated with turning points of the Alfven continuum, 
since these tend to dominate the oscillatory properties when one computes the time 
evolution of systems with continua (Levin, 2007; Sotani et ah, 2008). 

In fact, for a star with a magnetar strength held, crustal vibrations and core 
vibrations should couple together on very short timescales (Levin, 2006, 2007). 
Considering them in isolation, as described above, is therefore not appropriate. The 
current viewpoint, based on more detailed modelling that takes into account the 
magnetic coupling between crust and core, is that the QPOs are in fact associated 
with global magneto-elastic axial (torsional) oscillations of the star (Glampedakis et al., 
2006; Lee, 2008; Andersson et al., 2009; Steiner & Watts, 2009; van Hoven & Levin, 
2011, 2012; Colaiuda & Kokkotas, 2011, 2012; Gabler et al., 2012, 2013; Passamonti 
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& Lander, 2013, 2014; Asai & Lee, 2014; Glampedakis & Jones, 2014). However since 
magneto-elastic oscillations depend on the same physics described above, albeit now in 
a coupled system, they have frequencies in the same broad range as the simple estimates 
given above. Current magneto-elastic torsional oscillation models can thus in principle 
explain the presence of oscillations at frequencies of 155 Hz and below. 

Until very recently, however, it appeared that there was a signihcant problem with 
the higher frequency QPOs. This is because although there are crust shear modes 
in this frequency range, they should overlap with the various Alfven continua (there 
are no gaps between the harmonics of the continua as the frequency increases). As a 
result, the coupled oscillation should damp very rapidly, on timescales of less than a 
second (van Hoven & Levin, 2012; Gabler et ah, 2012). The data analysis, however, 
indicated that the oscillations persisted for up to ~ 100s (Watts & Strohmayer, 2006; 
Strohmayer & Watts, 2006). Various solutions to this problem have been explored, 
including coupling to polar modes (Lander et ah, 2010; Lander & Jones, 2011a; Colaiuda 
& Kokkotas, 2012), and resonances between crust and core that might develop as a result 
of superfluid effects (Gabler et ah, 2013; Passamonti & Lander, 2014). It is clear from 
these studies that superfluidity in particular can have a large effect on the characteristics 
of the mode spectrum: and since superfluidity is certainly present in neutron stars, mode 
models must start to take this into account before we can make hrm mode identihcations 
(Fig. 24). 

However the debate over this issue also exposed the fact that the initial data analysis 
did not actually test whether the signal could also be there in much shorter data 
segments, more consistent with the theoretical predictions. Huppenkothen, Watts & 
Levin (2014) have since re-analysed data for the 625 Hz QPO in the SGR 1806-20 giant 
flare and found that the data are in fact consistent with a short-duration signal that 
damps and is re-excited several times (rather than a long-lasting low-amplitude QPO). 
What might cause late time excitation and re-excitation remains an open question, and 
is relevant to the lower frequency QPOs as well since several seem to appear only late 
in the tails of the giant flares. Aftershocks may play an important role in exciting and 
re-exciting the QPOs that we see, and there may also be intrinsic delays in the process 
whereby vibrations are excited by the flare due to impedance mismatching between the 
different components of the star (Link, 2014). 

Since giant flares are very rare, there have also been efforts to search for seismic 
vibrations in the much more frequent lower energy bursts:};. As discussed above, it is 
not yet entirely clear whether these bursts are caused by the same mechanism as the 
giant flares. However if they are, it is quite possible that they might excite seismic 
vibrations at frequencies similar to those seen in the giant flares, particularly if we 

f There have also been a number of searches for gravitational waves associated with magnetar flares 
and any associated starquakes or global seismic oscillations (Abbott et ah, 2007, 2008, 2009; Abadie et 
ah, 2011). So far only upper limits have been reported, but new analysis techniques are being developed 
for the next generation of detectors (Murphy et ah, 2013). 
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Figure 24 Simulations of QPOs of a magnetized neutron star with a solid crust and 
superfluid core, in General Relativity, from Gabler et al. (2013). The upper panels 
show the Fourier amplitude and the lower panels the phase. Both frequency and mode 
structure change as the held strength varies. The color scale ranges from white-blue 
(minimum) to orange-red (maximum) in the top panel, and from 6 = tt/2 (blue) to 
9 = —ti/2 (orange- red), respectively. The crust is indicated by the dashed black line, 
and magnetic held lines by the solid magenta lines (from Gabler et al., 2013, ©2013 
American Physical Society, reproduced with permission). 


are genuinely seeing global modes of vibration of the star§. Searching for QPOs in 
the smaller bursts is however complicated by the short, transient nature of the burst 
lightcurves themselves, and this has required the development of specially tailored 
statistical methods (Huppenkothen et ah, 2013; Huppenkothen et ah, 2014b). 

So far these techniques have been applied to several data sets. A sample of 27 
bursts from the magnetar SGR 0501-1-4516, using Fermi GBM data, made one candidate 
detection, but its signihcance was weak (Huppenkothen et al., 2013). A search of a larger 
sample of 286 Fermi GBM bursts from SGR J1550-5418, however, found signihcant 
QPOs at 93 Hz and 127 Hz after averaging together multiple bursts from highly active 
episodes (Huppenkothen et al., 2014). Similar analysis using RXTE data from the most 
burst-active magnetars SGR 1806-20 and SGR 1900-1-14 (the two sources for which 

§ The search for global seismic vibrations in small and intermediate magnetar flares is a core science 
driver for future hard X-ray and gamma-ray missions such as the proposed Large Observatory for X-ray 
Timing (Feroci et al., 2014). 
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QPOs have been observed in the giant flares||) led to the detection of a QPO at 57 
Hz after averaging together multiple bursts from SGR 1806-20 Huppenkothen et ah 
(2014b). These frequencies are in the range found in the giant flares, and the QPO 
widths are also comparable. It therefore seems plausible that they are instances of 
the same phenomenon. If these frequencies do indeed represent global magneto-elastic 
oscillations the implication is that such vibrations are excited not only by giant flares, 
but also by trains of shorter bursts. This is important information when we start to 
consider how modes are excited by the trigger mechanism. 

The analysis of SGR J1550-5418, however, also revealed a QPO in a single burst, 
at a much higher frequency of 260 Hz. In addition to being in a different frequency 
band, this QPO was much broader than those seen in the giant flares and had very 
high fractional amplitude. If this is a magneto-elastic oscillation mode, then it is of 
interest since models predict that modes in this frequency range should die out on 
timescales comparable to the duration of short bursts. This could explain the observed 
low coherence, since broad width is a natural consequence of a rapidly exponentially 
decaying signal. This signal could, however be something quite different, such as a 
plasma instability associated with magnetic reconnection (Kliem, Karlicky & Benz, 
2000) or a local oscillation in a smaller, temporarily decoupled, cavity (Huppenkothen 
et ah, 2014). In this case it may be a hngerprint of the burst trigger process. Variability 
in the impulsive phase of the giant flares has previously been suggested, but dead time 
and saturation effects strongly distort timing analysis for the very brightest events 
(Barat et ah, 1983; Terasawa et ah, 2006). 

Another open question is how magneto-elastic oscillations couple to the 
magnetosphere and hence modulate the emission from the star. An important concern is 
that the fractional amplitude of the QPOs is in some cases quite high, and certainly much 
higher than the likely amplitude of any oscillations of the neutron star’s crust. Emission 
in the tails of the giant flares is dominated by radiation leaking from the trapped pair- 
plasma fireball. This emission is strongly beamed (giving rise to the strong rotational 
pulse), and thus in principle could act to amplify small surface vibrations, however 
analysis of the beams from the giant flares indicates that although the effect is real it 
is unlikely to be strong enough to explain the highest observed fractional amplitudes 
(D’Angelo & Watts, 2012). This suggests that there is some additional effect modulating 
the intensity of the emission: something that takes on added importance in the light of 
QPOs detected in the smaller bursts, where it is not clear that hreballs even form. A 
likely mechanism is a modulation of the optical depth to Resonant Gyclotron Scattering 
(see Sec. 3.3.1), via changes in particle number density (Timokhin, Eichler & Lyubarsky, 
2008) and/or magnetic held geometry (Gabler et ah, 2014). However the details of this 
process, and the interaction with the hreball, remain to be worked out fully. 

II Very few bursts from these sources have been observed in the period since Fermi GBM has been 
flying, and no giant flares have yet been observed in the Fermi era. 
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There is now a general agreement that the key observational phenomena that make 
soft gamma repeaters and anomalous X-ray pulsars so unique are well explained by 
the presence of a magnetar, a neutron star with ultra-strong magnetic held. In the 
absence of accretion from a binary companions and of large enough rotational energy 
losses, magnetic energy appears the only reasonable option to power both the persistent 
and bursting emission at the observed levels. Recently, thanks to the discovery of the 
so-called “low-B” sources, it has become increasingly evident that to make an “active” 
magnetar what matters is not (or not only) a large (> 10^^ G) dipole held, but a 
strong, residual poloidal component of the internal magnetic held. In this section we 
summarize the status of theoretical modelling, within the magnetar scenario, in the 
attempt to highlight what are, in our opinion, the issues which are basically settled and 
those which are still open. 

A clear picture has now emerged of where magnetars stand in relation to the 
population of isolated neutron stars, at least in terms of their observational properties. 
However, borders among the diherent classes are somehow blurry and the possible 
(evolutionary ?) links remain still to be fully understood. The existence of “low-B” 
and transient magnetars, implying that the number of highly magnetized neutron stars 
may be much larger then previously thought, makes it apparent that there may be a 
“birthrate” problem, unless there is an overlap between classes, meaning that objects 
that we see as observationally diverse are just neutron stars at different evolutionary 
stages, or magnetars can form through channels different from standard supernova 
events. 

A better understanding of magnetar’s formation path may help to bring clarity here. 
One of the biggest open questions is how magnetars acquire their super-strong magnetic 
held, as compared to those of the, apparently much more abundant, radio pulsars. 
Despite much effort, no dehnite conclusion has been reached as yet. Nonetheless, the 
very recent discovery that the magnetar in the Westerlund 1 cluster may have originated 
in a binary system points quite strongly towards a particular formation route involving 
massive binaries. Much theoretical attention has also been given to the idea that 
newborn magnetars are the central engines for gamma-ray bursts and newly developed 
models have had a great deal of success in explaining the “plateau” phase of the observed 
lightcurve. A conclusive proof for this seems now to rest with the observers, and may 
come in the next few years if a GRB is dehnitely linked to the gravitational wave 
detection of a compact binary interaction. 

A topic of the greatest relevance is the magnetic held conhguration in newborn 
magnetars. In particular how the held structure varies across core, crust, and 
magnetosphere; the balance between toroidal and poloidal components; and the small- 
scale structure of the external held. The state of the magnetic held at birth, and its 
subsequent evolution, are critical input for many aspects of magnetar physics. This 
is an active area of study that has also been given new impetus by the discovery of 
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magnetars with low surface dipole field strengths. There has been a good deal of progress 
in this area in recent years, although the various physical processes likely to affect the 
evolution of the magnetic held in the crust, where coupling of magnetic and thermal 
effects becomes important, remain challenging to model. 

There has been great progress in understanding the emission processes of magnetars. 
The existence of a twisted magnetosphere seems now to be generally accepted. Resonant 
cyclotron scattering onto pairs howing in a twisted magnetophere provides spectra 
which are in good quantitative agreement with the soft X-ray data, and also a natural 
explanation for the emission in the other wavebands (hard X-ray, optical/IR). However, 
a complete solution of the non-linear charge acceleration problem, including the various 
QED effects leading to photon splitting, positronium dissociation, and the fact that 
the twisted magnetosphere is expected to be strongly dynamic, has not yet been found. 
Another of the major open issues is the lack of a credible model of the interplay between 
crustal, atmospheric and magnetospheric emission, capable of explaining the broadband 
spectral energy distribution of magnetar sources. 

For the bursts, some aspects do seem clear. It is widely accepted, in the absence 
of a better model, that reconnection is likely to be required to explain the observed 
gamma-ray emission. For the giant flares, plasma ejection must take place to explain 
the radio afterglows, and a magnetically trapped pair plasma fireball seems the only 
viable hypothesis for the pulsed tails of the giant flares. The trigger mechanism for 
all bursts, however, remains unknown, as do the emission processes in the smaller and 
intermediate bursts. It will also be important to ensure that the bursts and persistent 
emission are considered as a whole: constraints on the magnetospheric structure and 
radiative transfer environment obtained from study of the persistent emission should 
be applied consistently, for example, when considering the radiation propagating in the 
aftermath of a burst. The detection of quasi-periodic oscillations in bursts has also 
provided new insight. The idea that these are caused by global seismic vibrations, 
excited by the burst process, is certainly the most plausible model put forward to date, 
although details of how the stars oscillate remain to be worked out. Far more work is 
required, however, is to examine self-consistently the conditions for excitation, decay, 
and modulation of the emission. 

Some theoretical predictions must await more advanced observational capabilities 
in order to be tested fully. NuStar (and possibly in the future ASTRO-H) is now 
offering the hrst opportunity to provide simultaneous data on the hard X-rays/soft X- 
rays turn over (at few tens of keV). ATHENA (the second large mission that will be 
developed by FSA, with launch in 2028) is the most important X-ray mission on the 
horizon and will have an unprecedented capability for compact objects and collapsar 
physics. A large area X-ray timing mission, such as the Large Observatory for X-ray 
Timing {LOFT, studied as a candidate for an ESA M3 mission), would also enable 
a large increase in capability to detect seismic vibrations in magnetar bursts and 
resolve both any mode splitting and their evolution on short timescales, two things 
that would both help to distinguish theoretical models. This may open the possibility 
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of performing, systematically, asteroseismology studies in neutron stars. The role of 
gravitational wave observatories in pinning down the mechanism behind short gamma- 
ray bursts, and the possible role of millisecond magnetars in that process, has already 
been described above. In fact, magnetars are wonderful candidates for detection by 
ground-based, long-baseline, interferometric gravitational wave detectors such as LIGO 
and Virgo. Magnetars are also powerful probes of the Galactic structure and the 
interstellar medium: pulsar timing arrays are also starting to be used to hunt background 
stochastic gravitational waves. 

Magnetar radiation is expected to be strongly polarized, and the polarization 
observables may also probe the so-called “vacuum polarization” effect, which is predicted 
by nonlinear QED, but has not yet been verihed experimentally. Future X-ray 
polarimetry experiments, currently under consideration for several small and medium 
missions (e.g. IXPE, a NASA SMEX candidate, and XIPE, an ESA M4 candidate) 
may therefore open a completely new window on our understanding of the radiation 
processes around magnetars and on the physics of matter and radiation in superstrong 
helds. 
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